Fonctions de Mathieu associées périodiques, équation différentielle de Mathieu associées 


solutions périodiques 


Je vais présenter dans ce chapitre une extension des fonctions de Mathieu, dont il est en général 
peu fait mention, tout du moins dans des articles récents. Ce sont les fonctions de Mathieu 


associées dont l'équation différentielle est la suivante : g''(¢)+ fa -2q Cos(21) er D et -o0 
Sin” (t 


D'autres fonctions de Mathieu associées intermédiaires sont liées aux deux équations 
u''(t)+ 2A Cotan(t) u'(t)+ (a — N —-2q Cos(2t))u(r) = 0 


différentielles équivalente : 
u''(t)— 2A Tan(t) u'(t)+ (a — N -2q Cos(2r))u(r) = 0 


Voici à ce jour de la rédaction de ce document les références que j'ai pu trouver dans la littérature 
scientifique à leur propos : 

- une article de 1921 de PHumbert : « On-Mathieu-Functions-of-Higher-Order » 

- un article de 1923 de E.L.Ince, « Associated Mathieu Functions », Proc. Edinb. Math. Soc., XLI 

- un article de 1926 de PHumbert, « XVII.—Some Hyperspace Harmonic Analysis Problems 
introducing Extensions of Mathieu's Equation » 

- un livret de 1926 de PHumbert, « Fonctions de Lamé et fonctions de Mathieu », pages 39 à 40 

- un livret de 1936 de P.Humbert, «Potentiels et prepotentiels », avec des mentions formelles mais 
non étudiées spécifiquement en pages 32, 33, 50, 51 

- 4 articles de R.Campbell, 1949, « Sur une expression remarquable des solutions de période 2kn de 
l'équation de Mathieu associée », 1950, « COMPORTEMENT DES FONCTIONS DE MATHIEU 
ASSOCIÉES POUR LES GRANDES VALEURS DES PARAMÈTRES », 1950, « CONTRIBUTION A L'ÉTUDE 
DES SOLUTIONS DE L'ÉQUATION DE MATHIEU ASSOCIÉE » et 1950, « Équations intégrales des 
fonctions de Mathieu associées et applications » 


Ces deux équations sont des équations de Hill soit des équations différentielles dont les coefficients 
du terme de dérivée première et de dérivée d'ordre O sont des fonctions périodiques de t. Ici la 
période est 2n, il est donc naturel de rechercher d'abord des solutions de ces équations de périodes 
JT OÙ 271. 


Ces équations différentielles sont intimement liées à l'équation des ondes sphéroïdales. C'est 
d'ailleurs une des voies originales pour leur introduction, tout comme pour l'équation de Mathieu 


qui en est aussi un cas particulier. 


Lorsque le paramètre A=0 ou A=1, l'équation de Mathieu associée devient l'équation de Mathieu. 


as Ç _ 24 Cos(2r) As) asd 


Lorsque le paramètre A=0, les deux équations intermédiaires de Mathieu associée deviennent 
respectivement des équations de Mathieu : À = 0" "()+(a—24 Cos(2r))(6)= 0 

u, "' (t)+ (a —2q Cos(2t)) u, (£) =0 
Voyons comment E.L.Ince dans son article véritablement introductif parvient à les obtenir à partir 
de l'équation des ondes sphéroïdales. 


Lien des fonctions de Mathieu associées (intermédiaires ou non) avec les fonctions sphéroïdales 


Une des formes de Liouville de l'équation des ondes sphéroïdales est obtenue avec le changement 
de variable x=Cos(t) : 


(1 T T a Ç -y?(1-x?) 1 hp x = Cos(t) 


o = Sin” (0170) > w''(t)+ (2u + 1) Cotan(t) w' (t)+ lo = u(u + 1)+ y’ Sin? 0) w(t) =0 
y(x) = Sin“ (0170) — w''(t)+ (1 — 2u) Cotan(t) w'(t)+ lo + ul = u)+ y? Sin°(t)) w(t) =0 


Une forme équivalente est obtenue avec y=Sin(t) : 
2 2 u 
G a 2x a + Ç H y0 x?) i z | y(x) =0 x= Sin(t) y(x) = (Cos? (e) 2 w(t) 


eN ei (1 + 2u) Tan(t) w'(t)+ lo — u(u + 1)+ y? Cos? (1) w(t) =0 
w" (t)- (1 — 2u) Tan(t) w'(t)+ lo + ul z u)+ y°Cos (£) w(t) =0 


Selon un article de 1923 de E.L.Ince, « Associated Mathieu Functions », le point de départ est la 
forme de Liouville x=Cos(t) suivante : 


1 1 
1 2 2 H+ Sur 
x=Cos(t)  p(x)=(sin ()) 4e) g'()+ w+ +11 Cos(2r) | ale 3) g(t)=0 
Posons A=p+> q= r a = © + Lr => e"O+fa+2a Cos(2t) AA Jet) = 0 


E.L.Ince réalise deux transformations supplémentaires : 


A 2A-1 
2 


a (in (DE n (= y) = (SOE e) (si) + 20 (0 = (sin (D) ale) 
= (SOE 2,02 va) (sr D) e) (sr) + 2,0) = (sir D)? ul) 


Ces transformations nous ramènent exactement à la première forme de Liouville proposée avec 
x=Cos(t), à savoir des fonctions solutions des équations différentielles suivantes : 


je "(z)+ 2A Cotan(t)u,'(t)+ (a -A +2q Cos(2t))u, (t) =0 
u» "()+ 2(1- A) Cotan(t) u,'(t)+ (a -(1— A) +2q Cos(2r))u, (t) =0 


Si l'on pose À =1_ u» le système d'équations est identique car l'équation est invariante par le 
2 


changement À 1-A< u >-+ On arrive donc à ces deux formes proposées par E.L.Ince : 
2AEL u 
ET re E nO One 
2 2A 


I)E (S E 2,0) = (sir OP 220) 
E 2A Cotan(t) u '(t)+ (a -4 +2q Cos(2t))u, (£) =0 
u,''(t)+ 2(1 = A) Cotan(t) u,'(t)+ (a = (1 = AY +2q Cos(2t))u, (£) =0 


Remarque importante : dans ce qui suit c'est le choix _ 1 _ u qui est retenu. 
2 


Remarque : sur cette généralisation de l'équation de Mathieu on remarquera que le paramètre q 
est de signe inverse dans l'article de E.L.Ince. En effet communément l'équation de Mathieu se 


note : y (t)+(a—2q Cos(2t))»()=0, et par extension l'équation intermédiaire de Mathieu 


associée : »‘'()+2A Cotan(t) y'(t)+ (a — À -2q Cos(2r)) »(1) = 0, et l'équation de Mathieu associée : 


gO a24 coste- AA D} 0 


Donc en revenant à la notation commune, il vient : 


Rae: ia à 
D) = (Sr (0) + 2, (0) = (Si OO) 2, (0) 
fe “(t)+ 2A Cotan(t)u,'(t)+ (a - À -2q Cos(2t))u, (t)=0 
u," (t)+ 2(1— A) Cotan(t)u,'(t)+ (a -(1- A} -2q Cos(24))u,(t)=0 
a (sin)? mO) v) = (sin O) + 200 = (sin (0) + E) = (sir D)? nle) 
g(t) = (sr (E)E u) > va) = (Sin (D) + 0) = (sin (D) + 0, ()= (Sin (D uzl) 


1- 1-2A 
4 


Construction des solutions périodiques des équations intermédiaires des fonctions de Mathieu 


associées forme x=Cos(t) 


Partons maintenant de ces deux équations proposées par E.L.Ince : 
u "(z)+ 2A Cotan(t) u'(t)+ (a — N +2q Cos(2r)}u, (£) =0 
u," (t)+2(1- A) Cotan(t)u,'(t)+ (a —(1- AP + 2q Cos(24))u,(t)=0 
Et recherchons la solution sous la forme de fonctions périodiques ou anti-périodiques de période r 


par des développements identiques en la forme à ceux employés pour les fonctions de Mathieu. 
Prenons par exemple une solution paire et anti-périodique en t de période n : 


1=+0 
u,()= D Ain Cos((21+1)t) 
1=0 


Injectons cette forme dans l'équation différentielle : 
Sin(t) (u, "(z)+ 2A Cotan(t) VA '(z)+ (a — N° +2q Cos(2r))u,(r))= (0) 


Aÿ.,Cos((21 +1) 


Donne sur le terme t) | 'expression suivante : 


AS aq Sinell -1)r)+(4+(1+21- AY -a)Sin(21 t)- (q +0 +21+ AŸ —a)Sin(2(+1)r)+ q Sin(2(+2)r)) 


Qui doit s'annuler une fois les termes de l'expression répartis sur les termes voisins du 
développement (ils sont tous en fonction Sinus dans l'équation différentielle de départ). Cela 
conduit à la relation de récurrence suivante à 4 termes: 


1=0>(a-(A+1})4^+(-a+q+(A-3F)4^ -q 4% =0 
1=1—>q A^ +(a-(A +3) -q) AN +(-a+(A-5F +q) A! +q 4^ =0 
l=2 >q A^ + (a (A +57 q)4ê+(a+(A-7F +q)4^ +q 4 =0 


q A4 +la-(21+1+AF -q) 44a +( a+(A—21-3) 4 q) A; —q AP =0 


, l= A 2 12 
Prenons un développement de la forme : u (t)=ce^ (t)= D ASCos(21 t), Par le même procédé, on 
1=0 
parvient à la relation de récurrence à 4 termes suivantes : 


1=0—(a-A-g)4 +(-a+g+(A-2Ÿ)A4" -g 4% 20 
1=1— q A + (a (A +2 -g) At +(-a+(A—4Ÿ +4) A4 +q 4% =0 
1=2- q Aù+(a-(A+4Ÿ -g) AS +(-a+(2-6) +4) 4! + 4! =0 


a A5 + (a-Gi+A) -a)4+ (bar (A 2-2) +4) 48, 4484 = 


Ce sont les deux récurrences que présente E.L.Ince dans son article de 1922 « Associated Mathieu 
Functions ». 


Cette article présente la récurrence sous une forme matricielle. Mais dans ce qui suit je vais 
délibérément choisir de modifier un peu cette présentation du système matriciel pour les raisons 


suivantes : 


- d'abord il convient d'inverser le signe du paramètre q afin d'être cohérent avec le standard 
d'écriture pour l'équation de Mathieu (au moins celui de NIST Handbook of Mathematical 


Functions, Mathieu Functions » 


- ensuite il est utile de faire apparaître plus clairement un problème de recherche de valeurs 


propres, et pour cela je vais changer de signe pour l'ensemble du système d'équations linéaires 


Alors la forme matricielle de la récurrence avec un développement de cosinus impair : 


(i+AŸ qg-G-A) -q 0 
q —g+(3+4A) g=-(5-A) -q = 
0 q -q+(5+A} q-0-47 -q 
0 q q-(21+1- AY 
Kas 0 q —q+(21+1+AŸ 
| 0 
E- CO T 
(SE Re A PA 
S(M-aB)A=0&(B'M-al}a=0 B'-|° | : i 
A 0 0 
00nn 0 


a 


0 0 
-1 0 
-1 

0 0 
i> j 


On voit donc qu'il s'agit d'un problème de détermination des valeurs et de vecteurs propres de la 
matrice B*M. Les coefficients du développement sont les vecteurs propres calculées à partir des 


valeurs propres déterminées par diagonalisation. 


Et la forme matricielle de la récurrence avec un développement de cosinus pair devient : 


-q+ 2q-(2-AŸ -q 0 
q -qg+(2+AŸ qg-(4-AŸ -q T 
0 q -qg+(4+A) q-(6-47 -q 
0 q Pe qg-(21-A) 
ii 0 q —q+(21+A) 
| 0 
[1 1 1 1 
lu 
S(M-aB)A=0&(B'M-ai)A=0 B'-|° ; : | 
is 0 0 
O O noon 0 


ai 


On voit donc qu'il s'agit également d'un problème de détermination des valeurs et de vecteurs 
propres de la matrice BM. Les coefficients du développement sont les vecteurs propres calculées à 


partir des valeurs propres déterminées par diagonalisation. 


Au 


Dans l'article E.L.Ince introduit les notations des fonctions associées de Mathieu sous la forme, que 
je change légèrement suivant le changement d'argument (préfixe « c » pour x=Cos(t) ou « s » pour 
x=Sin(t), correspondance avec les fonctions de Mathieu ce et se lorsque A=0, i lorsqu'il s'agit d'une 
équation intermédiaire de Mathieu associée, et l'inversion du signe du paramètre q->-q déjà 


évoquée) : 


UE —2q Cos(2t)- A 


A l=+%0 


g(t)= cce;(t,q)= (Sin? o) ccei!, (t,q)= (Sin? QA 


gli)= ecet (4) = (SOE cceit al, a)= (S0 


A (q) Cos(21 t) 


1=0 


Arala ) Cos((21 F 1)z ) 


=0 1=0 
I-A 1-A 1=+00 
glt)= esez (tq) = (Si (0) ? ccei," (1,4) (sin (O) X 45 (a) Cos((21+1)1) 
1=0 
1-A 1-A 1=400 


gle)= cseh, a (t,q)= (Sin? (0)? cceis,"(,4)= (sin (0)? X A*a) Cos(21 t) 


1=0 


Évidemment lorsque A=0, les développements sont effectivement d'une forme équivalente aux 4 
types de solutions périodiques de l'équation de Mathieu : 


A=0=> g'(t)+(a —2q Cos(21)) g(t)= 0> 


1=+0 
cce, (t,q)= ccei,,(t,g)= ce, (t,g)= X, 45(g) Cos(21 t) 
1=0 
l=+0 
cce, alt, q) = cceis, alt, q) = cent, q) = Anla ) Cos((21 az 1)t ) 
1=0 


1=+0 


ces 2(t q) E Sin(t) cce, lt, q) = se,,2(4) g)= = Sin(t) ox Arala ) Cos((21 +1) 1)z =È Bal )Sin((2 + 2) t) 
1=0 1=0 
csel (t,q)= Sin(t) ccei, (t,q)= sena (t,q)= Sin(t) OX) ) Cos((21) = 5 B,(q)Sin((21+1)t) 


1=0 


Construction des solutions périodiques des équations intermédiaires des fonctions de Mathieu 


associées forme x=Sin(t) 
Parallèlement étudions l'autre forme de Liouville obtenue avec x=Sin(t), à savoir : 
d’y(x) _, dy(x) p? i 
(1 x?) ne 2x de + œ©+y°(1 x?) TE y(x)=0 x = Sin(t) 


Au q= r a=w+ +7 
. = Cos (w(t) => w''(t)— 2A Tan(t) w'(t)+ (a — À — 2qCos(2t)) w(t) =0 


y(x) = Cos“ (rw) w'(t)- 21- A) Tan(t) w'(t)+ (a -(1- AŸ -24Cos(2t))w(t)=0 


On note qu'il suffit de réaliser le changement d'argument t->n/2-t pour obtenir la même chose à 
partir des équations précédentes. On peut donc facilement trouver la récurrence d'une solution de 


la forme „(t)= SA, Sin((21 +1)z): 


Partant du développement : 5 A Cos| (Z - Nen AA Sin((21+1)9, il suffit 


d'alterner les signes des coefficients de la matrice. Soit sous forme matricielle on obtient le système 
suivant : 


(+A? q+8-47 q 0 7 0 1 1 0 0 o1Y 4 
qd. (G+GFAŸ + GA’ q æ O 1 1 0 A, 
0 q g+(5S+AŸ g+(7-AŸ 0 | 10 0 1 1 0 |) 4 
0 q ae g+(21+1-A) 0 o Où OI se 
0 q g+Q@i+14A) E 0 0 1 a| 4u 
| o0 | {0 o0 oL 
[1 -1 1 -1 1] 
0O 1 -l 1 1 
e (M-aB)A=08(B7M-a1]jaz0 B° 9 1 oe eop] a … 
z O ON ES 7 18, =(-1)7 
RS D Le 
Où 406 ui o iia O 


Par une linéarisation trigonométrique, on arrive au même résultat. 


A .` 2 . l=+0 . 
De la même manière la récurrence d'une solution de la forme : „u(t)= S ACos(21 t} solution de 
= 


l'équation :#""()—2A Tan(s)u'(e)+ (a-~ —24Cos(2r)}u(r)= 0, est construite comme suit : 


[q+ 2q+(2-AŸ q 0 
q g+(2+A) g+(4-A) q T 
0 q g+(4+AŸ (66-47 0 
0 q g+(21-AŸ 
o 0 q (21+A) 
| 0 
[1 SR 1 
0 1 -=l 
&(M-aB)}A=0&(B'M-alla=0 B=’ ; : 
Fe. do MO 
o o 


0 [1 1 0 0 0]\ T4 
01 10.. A, 
_ 0 0 1 1 0 …|||4, 
— 4 = 0 
0 0 O … … 0 K 
a 0 0 A 
o 0 o …| 
-1 1] 
1 1 
—1 it b, =0 i>j 
1 =|b,.| = 1)” 
... EIEN 
O lon 
oon) 


Donc pour l'équation différentielle g" (t)+ (a — 2q Cos(2t) 


types de solutions (notation préfixe « s» pour x=Sin(t), 
Mathieu ce et se lorsque A=0) : 


gle)= sce (t,q)= (Cos”()}? sceif (1,4)= (Cost) 
gle)= sse, a(t, 4)= (Cos” (9) sceit a (.4)= (Cost) 
gle)= sse, (1.9) (Cos (D)? seei, (1.4)=(Cos°( 


g()=sce,. (,q)= (Co? (A) scei} (t,q)= (Cos?( 


A) g(t)=0; on peut construire 4 
Cos (£) 


correspondance avec les fonctions de 


A Al 


2 


+00 


À; (q) Cos(21 t) 


M 


S> I 


EROR CEDD 


1-A = 


DE > Arala 


1-A l=+%0 


D) Dar 


) Sin((21+1)£) 


1-A 


q) Cos (21 (21t) 


Là encore lorsque A=0, les développements sont effectivement d'une forme équivalente aux 4 types 


de solutions périodiques de l'équation de Mathieu : 


sce,, (t, q) = scei, (t q) = Ce (e, q) T 


A=0> g'()+(a —2q Cos(2t)) g(t)= 0> 


ssena lt, q) = ssein a (tq) = Sezal t ,q)= 


sce, l q) z Cos(t}scei,, (t q) = Cerny (e, q) = Cos(s)5 


l=+0 


> 


1=0 


4()Cos(21 1) 


S si (osile): 


1=0 


1=+0 


se (t, q) = Cos(t)ssei,,(t, q) = Seana (t, q) 5 Cos() A Arla )Sin((21 + 1)£ )= DB, (q)Sin((21 + 2) t) 


1=0 1=0 


B;(g)Cos((21 +1)1) 


l=+0 
4,(g)Cos(21 1)= 


1=0 


1=+0 


1=0 


Construction des solutions périodiques des équations des fonctions de Mathieu associées (forme 


x=Cos(t)) 


Revenons maintenant à l'équation différentielle originale : o" (t)+ (a — 2q Cos(2t) 


Il est possible d'obtenir les récurrences des développements précédents : 


eae CIE “ > 44 (q) Cos(21 ) 
cceh,. a(t,q)= (Sinte) Z 45 ala) Cos((21+1)) 
te [o, z] 1=0 


cseh.(.9)= (SM ^ FS B4 (q) Sin((21+2)1) 


esedalta)= (Sin) ^ D BA (a) Sin((21+1) 1) 


En injectant la forme directement dans l'équation différentielle : 


g" (t)+ (a —2q Cos(2r) 


MAD pe 


Sin?’ (t) 
Et dans ce cas on obtient des récurrences à 5 termes. 


Pour le développement en cosinus pairs (forme x=Cos({t)) : 
1 


(~ -q-a)4 + a+29-(2-A) 4, 4A, 0 


N 


1 
2 


(a+2q vase a+ A ajasta (4-4, TA, 0 


qa + la +29 (+4) )4,+(16+A -q a+ (a+ 24 (6-A} 4, 44 0 


1 1 1 1 
tat (a+2q-(21-2+ AP Mna +42 +A -q-a)4 + > (+29-(21+2-A) hna -5da =0 


2 2 J E 
oea 4 0 : ol Ji i 
2 2 2 Ja 2 
2 
2-À 43,4% 2q-(4-A) q j NE 
2 ne 2 2 
g 2-C+A) Lo i 24-(6-A) q o | |0 _1 
2 2 2 2 2 
q 2q-(21-2+AŸ PT 2q-(21+2-AŸ q 
2 2 2 2 
0 RE 


EIE H)}A=0&(H'M-a1)A=0 


A(A-1) 


Sin’ (t) 


Jo- 0 


Pour le développement en cosinus impairs (forme x=Cos(t)) : 


(ta aja +sla+a B Aÿ}4, q4 = 0 


le+a-(+A)}4 +0 q+A a+ la+20 (s-A) 4, 


1 
JIE 


a+ la+2g (3+A)°)4,+(25+ n° -q al +la+ 24 (7) 4, 


1 1 1 1 
z ors i slez (21 1+ AP a +O} +A -q au +5 +24 (21+3 AP Mns 7 dés © 0 
+A}  g-6-A} q i ape 4 
2 2 2a 2 2 
azra pin 2-64) q ý a EE 
2 l 2 2... 2 
ao -BHA ag 24-(7-A) _q 0 pre 
2 2 2 2 -a 2 
e.. e ... 3 er. .. 5 srt 
1 a A EVAN A 


1 
7% 


0 


&(M-aH)}A=0e(H'M-a1)A=0 


Pour le développement en sinus pairs (forme x=Cos(t)) : 


le AŬ u a+ (a+ AP ]8, 298, 0 


Slaa-(-aP)a is g+ A -4A alB,+(a+29 (6+Aÿ }B, TaB, 0 


Obs + (+20 (4 A}, +66+n 4A-q a)8,+(u+29 (A 598, =0 


l l l l 
EE (1-A¥ B, +21 +2} -g+ N-44 aBa +5024 (+444) BB =0 


22-A)Ÿ-9 24-(4+A) q j 
T J 2o 
2q e-s) 16 LA 4A 24-(6+A) 3 0 
2 2 
2 2 2 2 
q 2q-(21-A) (+0 -q+ A -4n 2q-(21+4+A) 
2 2 2 
0 


© (M-aH)B=06 (H'M-a1)B=0 


Pour le développement en sinus impairs (forme x=Cos(t)) : 


: bi A) -44 5a)B + (c+34 G+aÿ}B, TaB, 0 


ln +0 qg+A-4A a)8,+>(a+29 (af }8.-2u8, 0 


taB + (a+ 2 (G3-Aÿ )8,+(25-4+ 8-41 alB,+(e+24 (7+aÿ 8, “a 0 


1 l 3 I l 
Bas +5 la+29 (21-1 aÿ}B, +(21+1) q+A -4A alex + la +29 (1+3+A} Bas ab =0 


3(1-A) -44  3q-(G+A) q : à po 3 
Din 2 2 2 2 
AE Ge te EN d 0 mi SE 2e 50 
2 baf 2 ; 2 2 
q 2q-(G3-A SA 2q-(7+A) q 0 -2 1 -9 
; 25-q+A-4A = 0 2y À 7 
i 4 29-(21-1-AŸ - j -434A g : ee í zŠ 


&(M-aH)B-0(H'M-41)B-0 


J'ai vérifié que les divers systèmes linéaires obtenus même quand ils sont finis, sont strictement 
équivalents à ceux obtenus avec les récurrences à 4 termes définies précédemment. De même que 


pour A=0, on retrouve strictement les récurrences à trois termes des 4 types de développement des 
fonctions de Mathieu classiques. 


solutions 


+(a-24 coster- MAD 


Quant aux de l'équation originale de Mathieu associée 


ko la correspondance avec les fonctions de Mathieu est 


Sin? (£) 
également obtenue lorsque A=0, sous la forme matricielle qui est exactement celle des fonctions de 
Mathieu , notamment pour les ordres impairs : 


1 1 1 1 … 1 1 q4-9 -q 0 0 
O 1 1 1 … 1 q 9-q q-25 -q 
B-M - 0 0 1 1 0 q 25-q q-49 R Ka 
0 0 … 1 0 q g-(21+1) 0 
7 0 0 1 ni 0 ga  (21+1) -q 
|0 0 | 0 
[1+4 q 0 0 l+g-a q 0 0 oll 4 
q 9 q 0 ii q 9—a q 0 ne- À, 
SBM- 0 q 27 0 EE 0 q 27-a q 0 io A; 
0 q q 0 0 q q 0 sas 
” 0 q (21+1) 0 q (2l+1f-a | | Ara 
| 0 | 0 Les 


Et pour les ordres pairs : 


[1 1 1 1 1] [-q4 2q-4 -q 0 0 
O 1 1 1 1 q 4-q q-l16 -q 
NE 0 O0 1 1 0 q 16-q q-36 0 , 
0 O0 … 1 ma 0 q q—-4 0 
je 0 0 1 3 ce 0 q 4P -q 
|o 0 0 | 0 
[oO 2qg 0 0 0 [—a 2q 0 0 oll 4 
4 q O0 de q 4-a q 0 À, 
BM = 0 q 16 q 0 le 0 q 16—-a q 0 À, = 
… 0 4q q 0 aie 0 q y q 0 si 
o q (2) … Re o q (2If-a ...| |4, 


RS ue A a ea a lle 


Développement direct des solutions périodiques de l'équation de Mathieu associée (forme 
x=Sin(t)), en fonctions sinusoïdales 


On peut également rechercher des solutions périodiques de l'équation différentielle de Mathieu 


associée TONE + 2q Cos(2t) i à) g(r)= 0, en injectant directement des développements en 
Cos‘ (t 


fonctions sinusoïdales. Plus précisément sous les diverses formes suivantes, cela conduit à des 
récurrence à 5 termes : 


Développement en fonctions de Cosinus pairs (forme x=Sin(t)) 


g(t) = (Cos (e) > A,Cos(21 t) 


A, Aa + (8, = xr'a) A +(7, = xla) An + (5, — x'a) Aora + Er Asa = O 
£ CXI 
MT B SE xr =1 
B, = (4 (21-2+AŸ) ;=4P-g+A 6, = (a +k2g—(21+2-2%)) 
1 1 
1=0 > (yo — xla) A, +(6, — xé'a) A +84, =0 Yo =5(~-4) xo =5 


1=1— (8 — x'a) A +(7, - x°a)A, + (5, — xi'a)A, +£,4; = 0 Vi =4-Žq+ À 


En introduisant M une matrice penta-diagonale et C une matrice tri-diagonale, le système 
d'équations linéaires devient le suivant : 


Yo o & 0 0 Xo Xo 0 | EM 

B On à & 0 0 A PA Xi 0 0|| 4 

a, Pr n € 0 0z xng o o 4 

0 a; B y; & €, 0j-aj|... 0 Feu x pA 0 || 4 |=0 
z A, 
L 0 in lil 


(M-aC)A=0(C'M-a1)A=0 


On voit donc qu'il s'agit d'un problème de détermination des valeurs et de vecteurs propres de la 
matrice CM. Les coefficients du développement sont les vecteurs propres calculés à partir des 
valeurs propres déterminées par diagonalisation. 


Développement en fonctions de Cosinus impairs (forme x=Sin(t)) 


g()=(Cos? O)E D ACos((2 +1)r) 


A 3 +(8, in ara) A, E (7, m xla) As + (5, r= xi'a) A3 + EA2145 = O 


£ sr 
a =86 =-5 xr =x" ta Ar =1 


B, = (ea (-1-AŸ) n =(21+1) -q+ -44A 5 = (24-(21+3+A)) 


1=0> (r — xla), +(5, - yÿ'a) A, +£€oAs = 0 Yo -26-1 -4q) ôo =-464-6+47) xe Le 


1=1— (8, - x'a) À +, — xla) 4, +(6, — xř'a)A, +e,4, =0 PB; = (63-(4-1Y) 


En introduisant M une matrice penta-diagonale et B une matrice tri-diagonale, le système 
d'équations linéaires devient le suivant : 


Pr So a 0 a & ON Pat ot 0 Ya] 

pı 71 ô, E 0 … 0 x xe xe 0 A À, 

æa, Pa n Ô & 0 … 0 X L ge 0 … 0 À, 

0 a; BP, y; 0, €, 0|-a|.… 0 +. %5 x O0 … A |=0 
T Aan 
LO A L L 


(M-aC)}A=0(C'M-a1)A=0 


On voit donc qu'il s'agit d'un problème de détermination des valeurs et de vecteurs propres de la 
matrice C'M. Les coefficients du développement sont les vecteurs propres calculés à partir des 
valeurs propres déterminées par diagonalisation. 


Développement en fonctions de sinus pairs (forme x=Sin(t)) 


g(t)= (Cos? Oa F BaraSin((21 +2)ż) 


&, Bai + (8, = x'a) Ba + (7, z; xla) Ba + (5, - x'a) Bora + EBan =0 


3 E! 
= =-T xr x 2 x: =1 
1( 2) SaN KA LC 2) 
B =-3 2q—(21-—A) y, = 44 +1) -q +A -4A ner 2q-(21+4+A) 


1=0 > (yo — xla) B, +(6, — xia) B, +£B=0 y, =(A-2Ÿ + 


En introduisant M une matrice penta-diagonale et C une matrice tri-diagonale, le système 
d'équations linéaires devient le suivant : 


Vin Oy Ey Où ee ae 0 Ho A a au te ae TO B, 

AEn à € 0 … 0 o Do æ NO se Ol B 

a, Pa n Ô & 0 … O OA DE AD 5 Bs 

0a a a 6 Ol-a.. 0 g gR g ON B |=0 
Bana 

[0 ] Lil 


(M-aC)}B=-0«(C'M-a1)B=0 


On voit donc qu'il s'agit d'un problème de détermination des valeurs et de vecteurs propres de la 
matrice C M. Les coefficients du développement sont les vecteurs propres calculés à partir des 
valeurs propres déterminées par diagonalisation. 


Développement en fonctions de sinus impairs (forme x=Sin(t)) 


g(t) = (Cos? (e) > B, Sin((21+1)r) 


Br 3 +(8, = xa) Bari + (7, = xla) Bara +(5, a xřa) Bar3 + EBan = 0 
1 


a =e =-T EDR z; =1 


B, a (21-1-AŸ) n =(21+1) -4+4 5 =->(4-(21+3+A)) 


1=0 > (yo — zla) B, +(6,-xi'a)B, + £B, =0 vo =+} 5 =-żl-6-17) x => 


1=1—> (8, - x;'a)B, +i - ľa) B, +(6, — xi'a)B, + £,B, =0 B, = -(1+A}) 


En introduisant M une matrice penta-diagonale et B une matrice tri-diagonale, le système 
d'équations linéaires devient le suivant : 


BA do & 0 0 | Xo xo 0 | B, | 
B On à & 0 0 RS D x 0 0 B, 
3 Pa n Ô & 0 … 0 1: x L 0 0 B; 
0 a B y; & €, 0j-al... 0 ls PA pA 0 … B, |=0 
mao mia ME Hi vi dé E Bara 
OS ce NS es tn e a 2 sie] 


(M-aC)}B=-0< (CM 18 = 0 


On voit donc qu'il s'agit d'un problème de détermination des valeurs et de vecteurs propres de la 
matrice C'M. Les coefficients du développement sont les vecteurs propres calculés à partir des 
valeurs propres déterminées par diagonalisation. 


Impossibilité de la construction des solutions périodiques des équations des fonctions 


intermédiaires de Mathieu associées forme x=Cos(t) 


Voyons si l'on peut directement construire pour l'équation différentielle intermédiaire : 


u''(t)+ 2A Cotan(t) u'(t)+ (a — N —2q Cos(2t))u(t) =0 


des développements périodiques possibles à la manière des fonctions de Mathieu, sous la forme 
SE BA aSin(@1+1)¢) ou sous la forme >. BSin((21)r). Pour t=0, ces deux développements 
1=0 1=0 


s'annulent, ainsi bien évidemment que leurs dérivées secondes. Les dérivées premières sont a priori 
non nulles. Si nous injectons ce résultat dans l'équation différentielle alors : 


u'"'(t)+ 2A Cotan(t) u'(t)+ (a — N -2q Cos(2r)}u(s). =2A Cotan(t) u'(t ) =0 


t=0 


Comme A<>0 (sauf pour les fonctions de Mathieu) alors la solution ne peut qu'être identiquement 
nulle. 


De même avec l'équation différentielle “'()-2A Tan(t)u'(t)+(a -^ —2g Cos(2r)}u(e)= 0 par 

simple changement d'argument, on en déduit qu'il n'est pas possible de construire des solutions de 
1=+00 1=+00 y: . 

la forme S Aĝa Cos((21+1)r) OU de la forme S B4 Sin((21)r) excepté dans le cas des fonctions de 
1=0 1=0 


Mathieu, soit lorsque A=0. 


Résumé des 


rincipaux résultats 


our la construction des fonctions de Mathieu associées 


périodiques, avec des récurrences à 4 termes : 


ge e(a- 24 Cos(2i) 


© 
O © = m. 


O Om m. m. O O e m. m. O O me m. m. 


O O e m.e m. 


A 1) 


Snt }«o- 0 


glt) Sp (4,q)=( Sin e)“ > AS \(g) Cos (21 t) 
g(t)= ccenn (t,q)=( Sin O) DEN (4) Cos(( 21+ 1) t) 
glt) = esez, (t,q)= (Sinf) ^ D (4) Sin(r)Cos((21 +11) 


g(t) = cseza (t, q)= (Sine) ^ 3 47^ (q) Sin(t)Cos(21 t) 


[A+AŸ g-(3-AŸ -q 0 ni 0 
q —-qg+(3+A) q-(5-A) -q X 
0 q —g+(5+AŸ qg-(-A) -q ” 
0 q —(21+1-A) 0 
p 0 q —q+(21+1+AŸ 
0 
CHEN 2q-(2-AF -q 0 o] 
q —q+(2+AŸ q-(4-AŸ -q x 
0 q —q+(4+A) qg-(6-A) -q x 
0 q 7 -21-A ol 
ni 0 q —q+(21+A) 
0 = 
(2-AŸ g-(2+AŸ -q 0 : 0 
q -q+(4-AŸ qg-(4+4) -q a 
0 q -q +(6-Af q-(6+A) -q 5 
0 q p —(21+AŸ 0 
7 0 q -qg+(21+2-AŸ 
0 
[-q+0-Af 2q-(1+AŸ -q 0 
q -q+(3-Af g-G+AŸ -q Ba 
0 q -g+(s-AŸ g-(5+4) -q 
0 q —(21-1+A) 
0 q —qg+(21+1-A) 


‘l-ax|. 


-al |. 


Résumé des 


te [o, z] g"'(t)+ Ê -2q Cos(2t)- 


Z sani 
à. à 
ab y at 
2 2 
g eki 
2 
| 0 
E al f 
I 
Aj aL 
2 2 
o -lı 
2 
|o 
ju 
2 
SR ER. 
2 2 
o -l 1 
2 
| 0 
+ +: 
F I 
Aa a 
2 2 
pati 
2 
|o 


rincipaux résultats 


our la construction des fonctions de Mathieu associées 
périodiques, avec des récurrences à 5 termes : 


ce(s ER E ocot ) 


l=1 


l=+% 
AA -1 ccenalt,q)=(Sinle) 2 Anala) Cos((21+1)r) 
\ D ko=0> : 
in =+0 
cse,2(,9)=(Sin(9)" Z Biala) Sin((21+2)1) 
1=0 
cseh.(t,q)=(Sinlt PEnio Sin((21+1)t) 
41 [2 Lo ] 
p Ea uka q i 3 E 
2 2 2. À 
2 
0 Poe EN) A 0 0 4 
à OA a (6-4) A 
1 2q-(2+A 2q-(6-A 4 
l o 4 24 2 q a 
: 3 ; 16+ A -q - ; alu les 
1 i . à 24-(21-2+A) 7, Ci:2-AŸ à a 
2 2 2 2 AP +N -q 2 2 A2 
4 L 0 | 
= | 2 2 ] 
PTa aba) 0 i RL 
2 p 2 2 ; Á 
0 1 (a 9-q+N e 1 0 0 4 
l 24-(3+4) 24-(7-AŸ 4s 
Lo . q q 2 q q 
25+4 -q ; 0 or) … 
L'on a 2-(2-1+4) A (1+3-AŸ 4 + 
LE 21+1) + A | An 
2 2 2 2 ( + ) + q 7 2 2143 
SE, 20 | 
a+ | 2 2 ] 
; 22-A)-4  24-(4+A) q r i 
JA 2 N 
0 2-(-A) i6 q+A -4A 2-(6+A) 0 0 
1 24-(4-AŸ 24-(8+A) 
END a q q z 2 q 4 
3 , 36-q+A -4A ; - dll 
7 al s e.. 3 e.. ... 3 ... 
Sotaa q 2q-(21-A) 2 ; 2q-(21+4+A q 
21+2 À -4A 
2 : 2 2 Gai 2 2 
ce ED | 
P [a(i-AŸ-49 39-(G+4) q i y 
9 ue 2 2o, 
0 ED ga 28 Gin) T 0 0 
1 2q-(3-A) 2q-(7+AŸ 
l o £ q q 2 q _4 
: ; : 25-q+A -4A > 5 OEE 
“i PA 1 ... ... 5 … ... ; ... . 
FA aS q 2q-(21-1-A) 2 ; 2q-(21+3+A) q 
l 21+1 A -4A 
2 2 7 2 CO 2 2 
S ] 


Cas limite q->0 pour un paramètre A entier 


On sait que dans ce cas l'équation des ondes sphéroïdales se ramène à l'équation différentielle des 
fonctions de Legendre associées. La forme de Liouville permet d'appréhender les solutions lorsque 
u demi-entier : 


(1 Jet 2x 2A) H pean- p)-0 AH x = Cos(t) y(x)= P” (x)=(1-x Jule) = (Sine) u(t) 
P | 2v+1) à 

u (9+ 2A Count |! 7 ) -x Joe 2p+1 X Jeep 

u"(t)+(1-2u)Cotan(t)u'(t)+ (v + uv- u+1)ult)=0 


On a donc à notre disposition toutes les solutions de cette équation, sous la forme de fonctions de 
Legendre associées d'ordre demi-entier et son développement sous une forme finie 


l=p 


uļlt)= 3° 4? cos (v B zp+l +21 +1) ‘); Dans notre cas le paramètre v est également demi-entier, aussi 


1=0 


des développements du type u(t)= S BA. ,Cos((27 +1) t) OÙ u(t)= S'BACos(2T £) sont équivalents 


l=—œ0 [= 


level RER ion p+ 27 =v-2Pt a= 24v- p 
en posant 1 ol i 
ig V—-p-— z V—-p+— 
T=i+ 2 T =1+ 2 
à 2 


La récurrence des coefficients du développement impair „(t)= ` BA. ,Cos((27 X 1) t) donne : 
l 
2. 2 
T=0 [eR (1 p Jar [GE so) ae 
2: 2 
T =1 [en (G p |as [020 (5 | p) s? 


=—0 


4 


2 2 
fe (27 +1 pJ |. (e2 (27 : sep) |a S(p-1X2v-2p+21+1)B4,, =-(1+1Y2v+21+3)B., 


Et la récurrence pour un développement pair de la forme „(t)= 5 BACos(21 t) donne : 
l=—œ0 
2: 2 
roð [e p) ae -(E Earp )at 
4 4 
2: 2. 
Ti (22 (2 py |a (C2 Grp} ) ai 


4 


2 2 
(E er 2) }25 (eP (27 + 2+ pf | Bh < 2(p-I{2v-2p+21+1) Bå =-2(1 +1(2v +21+3) BÅ., 


Ce sont donc les mêmes récurrences que pour les fonctions associées de Legendre d'ordre demi- 
entier. 


D'un point de vue matriciel l'équation #"()+2A Cotan(t)u'(t)+ (a A Ju(e)=0 


pour un développement impair « cosinusoïdal » à ce système linéaire : 


111 1l[0-p? -6+př 0 0 
dr Cl i| o -pf -6+př 0 
0 0 1 0 0 (5- py 
0 0 1 0 0 —(21+1+ př 
sa 0 T 0 0 (27+1-p) 
[Oo 0 0 || o 
|@-pf-a -12p -20p  —-28p IA 
0 (3-p} -a -20p -28p À, 
al 0 0 (5-pf-a -28p |4 
0 0 Le zi Ta 
oa 0 0  (2l+1-pf-a || 4 
E S 23 la 
Les valeurs propres de ce système sont 4={1-PŸ.(-pŸ.(5-p).… 
Et pour un développement pair « cosinusoïdal » au système linéaire : 
Per ate ds Arr e0es) 0 0 
O 1 1 1 1110 (2-pf -(4+p) 0 m 
0 0 1 1 0 0 (4—-p} 0 0 
On 0.. 1 0 0 —(21+ př 
je 0007 D lle 0 0 (2-p} 
[Oo 0 0 || 0 
| p-a —8p —16p —24p Ia] 
0 (2-pf-a -116p  —-24p À 
2 
als 0 (4-pf-a -24p | 4 
0 0 K | A 
0 0 O (2l-pf-a || 4, 
| 0 E 


a= pp} (4-p) 
Les valeurs propres de ce système sont PRE RES RES 


Â=-P correspond 


0 À, 
À; 
-all G =0 
0 P 
Au 
Sal 
=0 
eril 
o] [4,1 
À, 
ne dde 
0 és 
À; 
IE 
=0 


Plus généralement lorsque q=0, on peut revenir à une valeur réelle de A par la substitution A & -p. 
Alors on a les systèmes matriciels suivants : 


- pour un développement impair « cosinusoïdal » : 


| (1+AŸ -a 12A 20A  28A . | À 
0 (3+A) -a 20A 28A "i ll 
0 0 (5+A) -a 28A a -I 4 |o 
0 0 z se les 
a à 0 0 (rH4A) a 4 
|0 E Fe n . 11 
a={i+A),(3+A),(5+A),.  u'"(r)+2A Cotan(t)u'(r)+(2p+1X2p+1+2A)u(s)=0 


Récurrence à deux termes Aya = i Z . = eN DAN R 21-1 
+1+p\l- p- 


Et pour un développement pair « cosinusoïdal » : 


A —a 8A 16A 24A ” …|| 4 

0 (2+A)-a 16A 24A a … || 4 
0 0 (4+A) -a 24A g || 4 |o 

0 0 i Pe iE 

. 0 O  (21+A)-a …|| A4, 

| 0 JL 


a= A ,(2 T A A 2A Cotan(t)u'(t)+2p(2p+2A)uļlt)=0 


Récurrence å deux termes À = C a r N 2er à 1) A 
PNP- 


Revenons aux solutions liées directement aux fonctions de Legendre d'ordre demi-entier. A chaque 
x 2 Sai 2 x y ~ 
fois les premières valeurs propres P et (1- p) correspondent à un développement à un seul terme, 
2 2 
les deuxièmes valeurs propres (-p) et G-p) à deux termes, les troisièmes valeurs propres 


(=p) et (5-P) à trois termes et ainsi de suite. Or on a vu dans la méthode de construction des 
fonctions de Legendre associées d'ordre et de degré demi-entier que quelque soit la valeur de 
p>0 : 

- pour v=-(2p+1)/2 et v=1-(2bp+1)/2, le développement ne contient qu'un seul terme non nul. Pour 
v=-(2p+1)/2 c'est un développement en cosinus pair et pour v=1-(2p+1})/2 c'est un développement 
en cosinus impair 

- pour v=2-(2p+1)/2 et v=3-(2p+1)/2, le développement contient deux termes non nuls. Pour v=2- 
(2p+1)/2 c'est un développement en cosinus pair et pour v=3-(2p+1})/2 c'est un développement en 
cosinus impair 

- pour v=4-(2p+1)/2 et v=5-(2p+1)/2, le développement ne contient qu'un terme. Pour v=4- 
(2p+1}/2 c'est un développement en cosinus pair et pour v=5-(2p+1})/2 c'est un développement en 
cosinus impair 

- et ainsi de suite 

- pour p=0, le développement n'a qu'un terme. 


Il est clair que cette séquence correspond à celle que nous trouvons ici. On peut donc en conclure 
l'identification des valeurs du degré v, suivant la succession des valeurs propres établies dans le cas 
des développements pair et impair selon la correspondance suivante : 


a=p >v= => p v=-u&v+u=0 développement cosinus pair 

a = (1 pÝ > v= : Pp v=l-u<&v+u=l1 développement cosinus impair 

a= (2 př > v= : P v=2-u<v+u=2 développement cosinus pair 
= = a= (3 pÝ > v= : P v=3-u<v+u=3 développement cosinus impair 

a = (4 př > v= : P v=4-u<v+u=4 développement cosinus pair 

a= (5 př > v= : P v=5S-u&v+u=sS développement cosinus impair 


Cette liste de valeur propre n'étant pas ordonnée de façon croissante, la correspondance en valeur 
croissante est la suivante : 


Valeur n°1 a=0—v+u=p développement cosinus parité p 


Val R2 +u=p-!l 
aleur n af L=p 


développement cosinus parité p—l1 
Valeur n°3 v+u=p+l 


Valeur n°2p -2 à v+u=tl 
a=(1-p} > 


développement cosinus impair 
eee Valeur n°2p -1 v+u=2p-1l fe K 
2 
Valeur n°2 p 5 v+u=0 , | A 
=p > développement cosinus pair 
Valeur n°2 p +1 v+u=2p 


Valeur n°2p+2 a= (p 1) >v+u=2p+l développement cosinus impair 


Valeur n°2p+3 a =(p 2) >v+u=2p+2 développement cosinus pair 


Orthogonalité des fonctions intermédiaires de Mathieu associées 
Écrivons différemment l'équation intermédiaire de Mathieu associées : 


u''(t)+ 2A Cotan(t) u'(t)+ (a- A -2q Cos(2r))u(r) = 0 


S Sin (A) Z L [sin UE (a- A 29 Cos(2r))u(t)=0 


Les propriétés de périodicité des deux fonctions construites sont : 


cceis, (t + x) = cceiÿ, (t) 


ein (— t)= = ccei, () 


1=+0 
ccei?, (£) = D A} Cos(21 t)= = 
1=0 


ccei;,' (0) = = ccei}, (5) = ccei,' (x) = 0 


cceih, (+ 7)=-cceif,. (9) 


À E -A 
CCel nyi (- t) CCE, (£) 


1=+0 
cceir,. (£) = D AS Cos((21 + 1) t) => o 
1=0 


cceif,.,"(0) = ccei,, (z) = cceii,'(x) = 0 


Il suffit donc d'étudier la fonction sur l'intervalle [0,n] pour en déduire sa forme partout. 


Soit deux solutions de valeurs caractéristique a différentes : 


snoi sin (D (a A 24 Cos(21)}u (e)=0 

Sin 9 E {sin Ol (a, à -24 Cos(2r)}u, (0 =0 
= si uO [sn Ou (DL sr POP) a (000 
= sr oii si fn D (POP (a -aau =0 


En intégrant sur l'intervalle [0,7] il vient : 


RATES 20)" “aa à Sin? Ou (ok, (1) = 0 


Lorsque A>0 alors il vient toujours a, +a, => Í dt Sin?^(t)u (tu, (t)=0- Lorsque A=0 alors 


c'est le cas des fonctions de Mathieu et nous avons l'orthogonalité complète sur les fonctions de 
Mathieu de même parité sur [0,7] et de parité quelconque sur [0,27] : 


Aan F Am => [a Un (07% (t) = 0 
0 
Ann À mn P fa Urns (us (e) =0 
0 
2x 
a, £a, > fa u, (tu, (£) =0 
0 


Supposons A>0, il vient pour les fonctions de Mathieu associées la relation d'orthogonalité 


suivante : a, Ż4„, > far cce^ (t)cce^(t)= 0: 
0 


Pour la normalisation des deux types de fonctions de Mathieu associées, 


Partant de la relation précédente d'orthogonalité des fonctions de Mathieu associées sur [0,27] : 
2x 
a #4, > [a cce^(t)cce^(t)=0 
0 


m 


On aurait pu calculer la norme sur [0,27] à partir de l'expression : fa ceed =2 fat t (ece (OF , mais 
0 


cette dernière conduit à des intégrales malaisées à calculer : 


a 2 
fat Sin" p A^ (g)Cos(21 ) | 
; >> far Sin" (t)Cos(21 t)Cos(2m t) 


0 


fat Sin* Zo )Cos(( (rno) 
0 


On peut donc s'affranchir de la relation d'orthogonalité des fonctions et par exemple étendre la 
règle de normalisation des fonctions de Mathieu (pour A=0), à savoir que la norme sur [0,27] du 
seul développement de Fourier est normalisé à 7 : 


l=+0 


f dt | D A4 (q)Cos(21 à) - = 7 = {4 (4) + 


1=0 


l=+0 


Sla- 


1=0 


l=+%0 


fa(Z a )\Cos((21+1)t D) =r => Ll Aiala = 


1=0 
2x 2x 


Ja (sm) (ces (a) = 


0 


o. 


0 0 


Et de même : 


Ja (Facon) -x = aa) + 


1=+0 


Ja (Fasen) -r> F= 
Ja (co) Geest = ar (cas) (sea) = 


0 


dt (Cos? O) (sses., (t, q)) 


o.m} 


dt (Sin? (O) (cc: (e, a) = [ai (sin cs t ,4) ) = fat t (Sin?( V (setul, a) = 
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Expression des dérivées premières des fonctions intermédiaires de Mathieu associées 


J'appelle les deux fonctions périodiques intermédiaires de Mathieu associées, les seuls 
développement en cosinus : 
u'"'(t)+ 2A Cotan(t)u'(t)+ (a — N —-2q Cos(2r))u(r) = 0 


1=+00 
u(t)= ccei/,(t,q)= © Ala) Cos(2l t) 
1=0 


Î=+00 
u(s)= cceiñ n (t,q)= >, A ,.(g) Cos((21+1)r) 
1=0 
Les dérivées premières et secondes sont alors triviales : 


1=+00 
Z { cceif, (t,q)}= -X 21 A4 (q) Sin(21 t) 
1=0 


2 1=+00 
< {ccei, (z, a)}= -J4 Alq) Cos(21t) 
1=0 
1=+00 


Sec. (.a)}= -X 1+1) 45, (a) SnI +09) 


d : 1=+00 
Sa (ccein ta) = -D (1 +1) 45, la) Cos((21 +11) 
1=0 


Développement des solutions périodiques de l'équation intermédiaire de Mathieu associée lorsque 
q proche de 0 pour un paramètre À quelconque 


L'approche que j'utilise pour établir les premiers termes d'un développement en puissance du 
paramètre q pour l'équation intermédiaire de Mathieu associée des valeurs caractéristiques a est 
purement algébrique, mais ne visent pas à établir une formule générique : 


u''(t)+ 2A Cotan(t) u'(t)+ (a — N —2q Cos(2r)}u(s) = 0 


C'est une approche tout à fait similaire et pragmatique à celle de McLachlan dans son ouvrage 
complet sur les fonctions de Mathieu. 


Prenons un développement pair et sa valeur caractéristique sous la forme suivante : 


T1=+00 
cceis(t,q)= >. A} (q) Cos(21 t) 
1=0 


Supposons maintenant que la valeur caractéristique et les coefficients se développent en série de 
puissance de q, connaissant par ailleurs la valeur caractéristique lorsque q =0 (liée à celle des 
fonctions de Legendre associées d'ordre demi-entier) : 


a, (a)=(n+AŸ +q af +g az +q° ai +. 
Aj(a)= 4° qA +q? A? +q? A +. 


En injectant ces formes directement dans l'équation : 


Sin(t) u" (t)+ 2A Cos(t) u'(t)+ Sin(t) (a -A —2q Cos(2t))u(r) = 0 


Et par linéarisation trigonométrique, il suffit de résoudre un système d'équations algébriques : 
- d'annulation de chacune des composantes de puissance de q 
- dans la composante de puissance d'annuler le terme trigonométrique linéarisé 


Il se trouve que le système d'équations algébriques est finalement linéaire en supposant un certain 
nombre de coefficients de départ comme des paramètres de départ et non comme variables à 
résoudre. Voici le code Mathematica permettant d'obtenir la solution pour le développement pair : 


Clear[a, a, AA, w, Ordre, I, jA,A, LaListeOrdreQ, ListSin, ListeSolve, ListeVariables]; 
ListeSolve={}; 
ListeVariables={}; 
Ordre=8; 
wlt_] := Sum[AA[2 1] Cosf2 It], {l, 0, Ordre}]; 
a=À[0]+ SumfAppendTofListe Variables, À[I]]; A[I] q^, {l, 1, Ordre}]; 
A[O]=(A)"2; 
Dof 
AA/2 1]= Expand[ Sum[AppendTofListeVariables, Af2 I, j]]; A[2 l, j] qi, {j, 0, Ordre}]] 

, {, 1, Ordre} 
J: 
Expr=TrigReduce[TrigExpand{Factor[Expandl Sin[t] (w^\[Prime]\[Prime])[t]+2 A Cos[t] Derivative[1][w][t]+Sin[t] (a-A^2-2 q Cos[2 t]) w[t]]]]]; 
LaListeOrdreQ=Take[CoefficientList[Expr, q], Ordre+1]; 
Do[ 

CArrays= 

Normal[CoefficientArrays{LaListeOrdreQ[[I+1]], Table[Sin[(2 1+1) t], {l, 0, Ordre}]]]; 

lf[Dimensions[CArrays] == {2}, 

ListSin=CArrays{[2]]; 
Scan[If[Not[# === 0], AppendTo[ListeSolve, # == 0]] &, ListSin] 

1 

, {l, 0, Ordre}]; 
ListeSolve=Simplify[ListeSolve]; 
leSolve=SolvefListeSolve, ListeVariables]; 
Print{Take[leSolve[[1]], Ordrel]; 
Print{Take[leSolve[[1]], Ordre] /. À -> 0]; 


Il suffit d'appliquer ce code pour les valeurs successives A[0]=(A)^2 , A[O]=(2+A)"2, A[0]=(4+4)^2,... 
pour obtenir à l'ordre 2, 3 ou 4 (suivant la complexité des expressions ordres supérieurs non 
reproduite ici), les valeurs caractéristiques suivantes : 


2A 1+2A A 2A(1+2A) > 
1+A T (AJGA) (AFC+HAG+A)" 
(1+2A-21-61A +464 +76A +20N) ,. 
4{i+A)(2+A)(3+A)X4+A) i 
2A(A-1) 45+111A +334 -19A -10A , 
344A+A TT (i+A)G+A)(G4+A) 
. 2A(-819-2697A -2676A? -1242M -259A +3A° +10AŚ) 
(i+A)(3+A)(4+A)5+A) E 
à 2  2A(A-1 3(-450-1005A—575A? -294 +3744 +6) , 
on se 2+A)3+4)(5+4A)(6+A) la | 
_ 6A(A-1(-37350-55065A — 290504? — 6950 4° —638A* +23 +6A°) ; 
| (2+A)3+A)(5+A)(6+AYX7+A) à 
2A(A-1)  19600+24465A +8015A° +53A°-267A*-26À , 
FANTA) (4+AX5+A)(7+A)(8+A) 
2A (A —1)-3586800- 5310095 -3012625A? — 828940 A° 1125541 —5567A +235A°+26A') ; 
(3+AX4+AX5+A)(7+A)(8+A)X0+A) oo 


n=0 > a (q)= À 


n=1 > aå (q)=(2+A) 


n=3— a (q)=(6+A) 


un 


Remarque : si l'on suppose avec Mathematica une valeur propre de la forme a=(2n+A)°, il se trouve 
A,0 

que les variables de départ RTE [0...7] pour le système d'équations algébriques sont elles- 

mêmes en nombre variable, dès lors je n'arrive pas à obtenir un système d'équation 

intrinsèquement cohérent pour être soluble tout du moins dans Mathematica. 


Voici le code Mathematica permettant d'obtenir la solution pour le développement impair, en 
supposant également les développements en puissance de q : 


u''(t)+ 2A Cotan(t) u'(t)+ (a -A +2q Cos(2t))u(t) =0 
a= azla) 


l=+0 
ccei,.\(,g)= X Aia (q)Cos((21+1)t) 


A 2 Al 2 _A, 3 A3 
a;,4(q)= (2n +1+A) +q Arna FA din FO din te. 
A __ yA.0 Al 2 4A,2 3 A3 
Ain (g) = Ara + Ara +49 Ai +4 Aara te 


Clear{a, a, AA, w, Ordre, |, jA, A, LaListeOrdreQ, ListSin, ListeSolve, ListeVariables]; 

ListeSolve={}; 

ListeVariables={}; 

Ordre=5; 

wit_] := Sum[AA[2 1+1] Cos[(2 +1) t], {l, 0, Ordre}]; 

a =À[0]+ Sum{AppendTo({ListeVariables,A[I]];A[1] q^, {l 1, Ordre}]; 

A[0]=(1+A)"2; 

Dof 
AA[2 I+1]=Expand[ Sum{[AppendTo{ListeVariables, A[2 1+1, j]]; A[2 I+1, j] ^j, Ü, 0, Ordre}]] 
, {, 1, Ordre} 


Expr=TrigReduce[TrigExpand{Factor[Expand[Sin[t] (w^\[Prime]\[Prime])[t]+2 A Cosft] Derivative[1][w][t]+Sin{t] (a-A72-2 q Cos[2 t]) wlt]]]]]; 
LaListeOrdreQ=Take{CoefficientList{Expr, q], Ordre+1]; 
Dof 
CArrays=Normal{CoefficientArrays{[LaListeOrdreQ{[[l+1]], Table[Sin[(2 1+2) t], {l, 0, Ordre}j]]; 
lf[Dimensions[CArrays] == {2}, 
ListSin=CArrays{[2]]; 
Scan[If[Not[# === 0], AppendTo[ListeSolve, # == 0]] &, ListSin] 


1 

, {, 0, Ordre}]; 
ListeSolve=Simplify[ListeSolve]; 
leSolve=SolvelListeSolve, Liste Variables]; 
Print{Take[leSolve[{[1]], Ordre]]; 
Print{Take[leSolvelf1]], Ordre] /. A -> 0]; 


Il suffit d'appliquer ce code pour les valeurs successives À[0]-(1+A)12,A[0]=(3+A)"2, A[0]=(5+4)^2,... 
pour obtenir à l'ordre 2, 3 ou 4 (suivant la complexité des expressions ordres supérieurs non 
reproduite ici), les valeurs caractéristiques suivantes : 
2(A-1) 3(1+2A) >,  6(A-1K1+24)  ; 
24A | G+AÏ2+A) Ÿ C+AÏ+AU4+A) ? 
15(1+2A)(12-121A-4247+124 +4At) , 
4(2+A)(3+A)(4+A)5+A) 
2A(A-1) 160-464A-108A° +434 +14 , 
C+Aja+4)° (2+A)(4+A)(5+4) 
_2(A—1{-7680-22976A —18272A° -62724 -862A* +2345 +14A°) , 
| (2+A)(4+A)(5+AYX6+A) d 
2A(A-1) 6048 -9648A - 40204? -874° +1814 +22A° , 
(a+ AX6+ A)? (3+AX4+4)(6+A)(7+A) f 
, 2A(A -1-387072 -519840A — 2374084? -478084 3603A +1094 224°) : 
(3+AX4+4A)(6+A)(7+A)8+A) 
2A(A-1) 3(-17280-17760A -4792 +9A°+123A* +104) , 
(6+AX8+A)? (5+AX6+A)(8+A)(9+A) 
GA(A—1)-4285440- 5319168 — 2548544? — 595392 A? — 68614A —2833A +11IA°+IOA) , 
+ q + 
(4+AX5+AX6+A)(8+A)(9+A)10+A) 


n=0 > a (g)=(1 H A} 


n=1 > a} (g)=(3+AŸ 


n=2 > a! (q)=(5 A} 


n=3 > a! (q)=(7+A) 


Remarque : si l'on suppose avec Mathematica une valeur propre de la forme a=(2n+1+A)°, il se 


A,0 
trouve que les variables de départ ia le [0-7] pour le système d'équations algébriques sont 


elles-mêmes en nombre variable, dès lors je n'arrive pas à obtenir un système d'équation 
intrinsèquement cohérent pour être soluble tout du moins dans Mathematica. 


Développement des solutions périodiques de l'équation de Mathieu associée lorsque q proche de 0 
pour un paramètre À quelconque 


C'est exactement la même approche que j'utilise pour la fonction associée de Mathieu liée à 
l'équation différentielle suivante : 


ARE -24 Cos(2r) a) aen 


Sin? (t) 
Pour les 4 types de développements en sinus et cosinus, pairs ou impairs : 
coes, (t,4)= (Sine) ŽS AAA) Cos(211) aa} (4) 
1=0 
cc, (1,4)= (Sin) E A5(g)Cos((21+1)9) a= ahala) 
1=0 
cseh.(1.4)= (Sin) ^ È Bho Sin((21+2)r) a =b},.,(a) 
1=0 


cseh (9) (Sin) X A5 (a) S(r) a=bh. (a) 


On suppose alors que les valeurs caractéristiques et les coefficients se développent en puissance de 
q: 


az (a)=(2n+A) +q ai +q? ai +q° azy +. 

az alq)=(2n+1+AF +q azh +g? ahata? ai +. 
banolq)=(2n+2-AF +q biia +q? bara +q? bia +.. 
banlag) =2n+1-AF +q ban +a’ bia +q° baa + 
AA lq)= 44° +q AA" + q°AË? + qéAËS +. 

Anala) = Axa + q Ah +P ANA +GP ANA + 

Ba alq) = Bx +q Baia +q Bah +q? Bii +. 
Baala) = B3 +q Bin +q? BXA +q° Baat- 


Voici le code Mathematica permettant d'obtenir la solution pour le développement pair en 
cosinus : 


Clear[q, a, AA, A, w, Ordre, |, j, À, A, LalListeOrdreQ, ListCos, ListeSolve, ListeVariables]; 
ListeSolve={}; 
ListeVariables={}; 
Ordre=7; 
wit_] := Sin{t]"A SumfAA/2 I] Cos[(2 1) tj, {l, 0, Ordre}];a=A[0j+Sum{AppendTo{Liste Variables, À[I]]; A[I] q^, {l, 1, Ordre)}]; 
A[O]=(A)"2; 
Dof 
AA[2 I]=Expand[ Sum[AppendTo[ListeVariables, A[2 |, j]]; 
AL I, j] q^, {, 0, Ordre}]] 
, {l, 1, Ordre} 
L 
Expr=TrigReduce[TrigExpand{Factor[Expandl Sin[t]^-A (Sin[t]^2 (w^\[Prime]\[Prime])[t]+(Sin[t]^2 (a-2 q Cos[2 t])-A (A-1) ) wft])]]]]; 
LaListeOrdreQ=Take[CoefficientList[Expr, q], Ordre+1]; 
Dof 
CArrays= 
Normal[CoefficientArrays{LaListeOrdreQ{[[l+1]], 
Table[Cos/2 It], {l, 1, Ordre}]]]; 
lf[Dimensions{CArrays] == {2}, 
ListCos=CArrays{[2]]; 
ListO= Simplify[LaListeOrdreQ{[[I+1]]- ListCos . Table[Cos[2 It], {l, 1, Ordre}]]; 
If[Not[ListO === 0], AppendTol{ListeSolve, ListO == 0] 


Scan[If[Not[# === 0], AppendTofListeSolve, # == 0]] &, ListCos] 


] 

, {l, 0, Ordre}]; 
ListeSolve=Simplify[ListeSolve]; 
leSolve=Solve[ListeSolve, ListeVariables]; 
Print[Factor[Take[leSolve[[1]], Ordre]]]; 
Print[Take[leSolve[[1]], Ordre] /. A -> 0]; 


Il suffit d'appliquer ce code pour les valeurs successives A[0]=(A)^2 , A[O]=(2+A)"2, A[0]=(4+4)^2,... 
pour obtenir à l'ordre 2, 3 ou 4 (suivant la complexité des expressions ordres supérieurs non 
reproduite ici), les valeurs caractéristiques suivantes : 


2A 1+2A 2A(1+2A 
ou GrA}@+ A) | Sn an ll 
(1+2A( -21-614 +464 +764 +20N) 4. 
4{1+A)(2+A)(3+A)4+A) ii 

2A(A-1)  45+111A+33A-I9A° -10A , 
344A+A TT (AFGAV GA) 
2 A(-819-2697A 2676A? -12421 -259A +3A°+10A°) 

(1+A)(3+A)(4+A)5+A) La 
2A(A-1) 3(-450-1005A - 5754? -294° +37At+6N) , 
G+AX+4)° (2+AX3+4)(5+A)(6+A) 
6A(A -1 -37350- 55065A — 29050 A7 — 6950A? — 638A* +23A°+6A°) , 

(2+A)3+A)(5+A)(6+A)X7+A) i 
2A(A-1)  19600+24465A +80154A? +534? -267A -26A , 
(S+AX7+A) (4+AX5+A)(7+A)(8+A) 

_2A(A—1)-3586800-5310095A -3012625A? — 8289404 —112554A —5567A° + 235A° +26A') g 


(3+AX4+4AX5+4)(7+A)(8+AY9+A) 


n=1 > aå (q)=(2+A) 


n=2 > ak (q)=(4+AŸ 


n=3- a(q)=(6+A) 


Voici le code Mathematica permettant d'obtenir la solution pour le développement impair en 
Cosinus : 


Clear[q, a, AA, A, w, Ordre, |, j, À, A, LaListeOrdreQ, ListCos, ListeSolve, ListeVariables]; 
ListeSolve={}; 
ListeVariables={}; 
Ordre=5; 
wlt_] := Sin[t]^A Sum[AA/2 1+1] Cos[(2 I+1) t], {l, 0, Ordre}]; 
a=À[0]+Sum{AppendTo{Liste Variables, [1]; A[I] q^l, {l, 1, Ordre}]; 
A[0]=(1+A)"2; 
Dof 
AA[2 1+1]=Expand[ Sum{[AppendTo{ListeVariables, A[2 1+1, j]]; A[2 I+1, j] ^j, Ü, 0, Ordre}]] 
, {, 1, Ordre} 


Expr=TrigReduce[TrigExpand{Factor[Expandl Sin[t]^-A (Sin[t]^2 (w^\[Prime]\[Prime])[t]+(Sin[t]^2 (a-2 q Cos[2 t])-A (A-1) ) wlt])]}]]; 
LaListeOrdreQ=Take{CoefficientList{Expr, q], Ordre+1]; 
Dof 
CArrays=Normal{CoefficientArrays{[LaListeOrdreQ{[[l+1]], Table[Cos{(2 1+1) t], {l, O0, Ordre}]]]; 
lf[Dimensions[CArrays] == {2}, 
ListCos=CArrays{[2]]; 
Scan[If[Not[# === 0], AppendTofListeSolve, # == 0]] &, ListCos] 


1 

, {, 0, Ordre}]; 
ListeSolve=Simplify[ListeSolve]; 
leSolve=SolvelListeSolve, Liste Variables]; 
Print{Factor[Take[leSolve[[1]], Ordre]]]; 
Print{Take[leSolve[f1]], Ordre] /. A -> 0]; 


Il suffit d'appliquer ce code pour les valeurs successives [0]-(1+A)12,A[0]=(3+A)"2, A[0]=(5+4)^2,... 
pour obtenir à l'ordre 2, 3 ou 4 (suivant la complexité des expressions ordres supérieurs non 
reproduite ici), les valeurs caractéristiques suivantes : 


2(A -1) 3(1+2A) >»,  6(A-IK1+24) ; 
27A 1 GARA) CAGANHA 
15(1+2A)(12-1214 -42A +124 +44) , 
4(2+A)(3+A)(4+A)5+A) 
2A(A-1) 160—464A —108A°+43A°+14A > 
C+AK4+A)7 (2+A)(4+A)(5+4) g 
(A —1)(=7680- 22976A 182724? -62724 -862A* + 23A° +14A°) A 
| (2+A)(4+4A)(5+AYX6+A) 5 
2A(A-1) 6048-9648A—4020A°-87A°+181A*+22K , 
G+AX6+A)7 (3+AX4+4A)(6+A)(7+A) i 
2A(A -1 -387072 -519840A — 2374084? - 47808A? —3603A +109A°+22A) , 
(3+AX4+4A)(6+A)(7+AÏ8+A) i 
2A(A-1) 3(-17280-17760A —4792A° +9A° +123A*+10A°) , 
(6+AY8+A)° (5+AX6+A)(8+A)(9+A) 
_ GA(A-1)- 4285440 5319168 — 2548544 A? — 5953924? — 68614A* —2833A +1 11A% +10A7) ; . 
(4+A)5+4AX6+A)(8+A)(9+A)X10+A) - 


n=0 > a" (q)=(1 A} 


n=l > ak (q)=(3+AŸ 


n=2 > a (q)=(5 H A) 


n=3 > a} (q)=(7 H AF 


Remarque : notons que pour la valeur A=-1/2, certains développements de valeurs caractéristiques 
se simplifient : 


A =X — j 

asla) TE 
2(A-1) 1 

A(a)=(1+A) =-+2 

a (g)=(1+A) un dt 


Ce sont effectivement les valeurs exactes des caractéristiques de l'équation de Mathieu associées. 
Toutefois ces valeurs ne sont pas nouvelles, car ce sont celles qui correspondent à la solution 
algébrique la plus simple de l'équation des ondes sphéroïdales. En effet : 


1 y? 1 y? 1 u=1 
A=—— = a=0@0+—4 =—+2q => 
E ‘ 4 4 2 4 4 æ=0 


les solutions de l'équation intermédiaire de Mathieu associée sont : 
u(t)= A Cos(2,/qCos(r)}+ B sin(2 /qcos(r)) 


Cela correspond à une solution algébrique de l'équation sphéroïdale bien connue : 


RE : A Cos(2/ax)+ 8 Sin (ax) _ À Cos(y x)+B Sin(y x) 


11- x? 1- x? 


On peut imaginer que les développements pour des valeurs de A=-3/2, -5/2, -7/2 correspondent 
effectivement aux valeurs analytiques des autres solutions algébriques de l'équation des ondes 
sphéroïdales pour u=2,3,4,... 


Voici le code Mathematica permettant d'obtenir la solution pour le développement pair en sinus: 


Clear[q, a, AA, A, w, Ordre, |, j, À, A, 
LaListeOrdreQ, ListCos, ListeSolve, ListeVariables]; 
ListeSolve={}; 
ListeVariables={}; 
Ordre=6; 
wit_] := Sin[t]^-A Sum[AA/2 1+2] Sin[(2 1+2) t], {l, 0, Ordre}]; 
a=À[0]+Sum{AppendTofListeVariables, A[1]]; A[!] q^, {l, 1, Ordre}]; 
A[O]=(2-A)"2; 
Dof 
AA[2 1+2]=Expand[ Sum[AppendTolListeVariables, A[2 1+2, j]]; 
A[2 1+2, j] q^j, {j, 0, Ordre)}]] 
, {l, 1, Ordre} 


Expr=TrigReduce[TrigExpand[Factor[Expandl Sin[t]^A (Sin[t]^2 (w^\[Prime]\[Prime])[t]+(Sin[t]^2 (a-2 q Cos[2 t])-A (A-1) ) w[t])]]]]; 
LaListeOrdreQ=Take[CoefficientList[Expr, q], Ordre+1]; 
Dof 
CArrays=Normal[CoefficientArrays{LaListeOrdreQ[[l+1]], Table[Sin[(2 1+2) t], {l, 0, Ordre}]]]; 
lf[Dimensions{CArrays] == {2}, 
ListSin=CArrays{[2]]; 
Scan[If[Not[# === 0], AppendTofListeSolve, # == 0]] &, ListSin] 


, {l, 0, Ordre}]; 
ListeSolve=Simplify[ListeSolve]; 
leSolve=Solve[ListeSolve, ListeVariables]; 
Print{Factor[Take[leSolve[{[1]}, Ordre]]]; 
Print{Takel[leSolve[[1]], Ordre] /. A -> 0]; 


Il suffit d'appliquer ce code pour les valeurs successives À[0]=(2-A)"2 , A[0]=(4-A)"2, A[0]=(6-A)"2,... 
pour obtenir à l'ordre 2, 3 ou 4 (suivant la complexité des expressions ordres supérieurs non 
reproduite ici), les valeurs caractéristiques suivantes : 


2A  3(2A-3) , G6A(2A -3) : 
ra {A-3J(A-4)7 (a -3ÿ(A-4XA 5)? 
15(2A-3(-135+153A+18A°-28A +4) , 
4(A—3) (A-4)(A-5XA -6) 
2A(A-1) _ —675+495A +1054? -994° +144 , 
a-a- araa * 
2 A(-56025+81585A - 418204? +9210A° 537 107A +1445) à 
(A-3)(A-5)(A-6)A-7) i 
2A(A-1)  19600-17115A+2975A° +8574 -2914f +22A° , 
(asa 77" (A -4A -5F(A-7F(A-8) 
2A(A —1)(=1195600+1151815A —401030A? + 60690A? -2728A* 2414" +22A°) , 
(A—-4{A-5)(A-7)(A-8XA-9) i 
2A(A—1)  3(89690-26775 A +39274? +6014? -173A +104) , 
(A-7%A-9)" (A-6{A-7P(A-9F(A-10) 
6A(A -112819870 -12490317A + 47728594? —895608A° + 807641 —1957A5 — 181A% +10A') ; 
(A-5{A-6XA-7F(A-9F(A -10A -11) 


n=0b;(g)=(2-A) 


n=1b}(g)=(4-A) 


n=2 be (g)=(6- A) 


n=3 bi (a)=(8-A) 


En fait il s'agit des mêmes valeurs propres que pour les développements impairs mais en réalisant 
la substitution A->1- A. 


Voici le code Mathematica permettant d'obtenir la solution pour le développement impair en sinus: 


Clear[q, a, AA,A, w, Ordre, |, j, À, A, LaListeOrdreQ, ListCos, ListeSolve, Liste Variables]; 
ListeSolve={}; 
ListeVariables={}; 
Ordre=7; 
wlt_] := Sin[t]^-A Sum[AA[2 1+1] Sin[(2 1+1) t], {l, 0, Ordre}]; 
a=À[0]+ SumfAppendToiListe Variables, A[1]]; A[!] q^l, {l, 1, Ordre}]; 
A[0]=(3-A)"2; 
Dof 
AA/2 l+1]= Expand[ SumfAppendTofListeVariables, A[2 +1, j]]; A[2 1+1, j] q^j, {j, O, Ordre}]] 
, {, 1, Ordre} 
Expr=TrigReduce[TrigExpand{Factor[Expandl Sin[t]^A (Sin[t]^2 (w^\[Prime]\[Prime])[t]+(Sin[t]^2 (a-2 q Cos{[2 t])-A (A-1) ) wf{t])]]]]; 
LaListeOrdreQ=Take[CoefficientList[Expr, q], Ordre+1]; 
Do[ 
CArrays=Normal[CoefficientArrays[LaListeOrdreQ[[l+1]], Table[Sin[(2 I+1) t], {l, 0, Ordre}]]]; 
If[Dimensions[CArrays] == {2} 
ListSin=CArrays[[2]]; 
Scan[If[Not[# === 0], AppendTof[ListeSolve, # == 0]] &, ListSin] 


1 

, {, 0, Ordre}]; 
ListeSolve=Simplify[ListeSolve]; 
leSolve=SolvelListeSolve, Liste Variables]; 
Print{Factor[Take[leSolve[[1]], Ordre]]]; 
Print{Take[leSolvelf1]], Ordre] /. A -> 0]; 


Il suffit d'appliquer ce code pour les valeurs successives À[0]=(1-A)2,A/0]-(3-A)"2, A[0]=(5-A)^2,... 
pour obtenir à l'ordre 2, 3 ou 4 (suivant la complexité des expressions ordres supérieurs non 
reproduite ici), les valeurs caractéristiques suivantes : 


2(A-1) 2A-3 2.  2(A-I{2A-3) 3, 
AI E ATA MW JU-MAT 
(2A -3)( 60-339A +3944? —156A° +204) , 
4(A-2)(A-3)(A-4YA-5) 
2A(A-1)  —160+80A+84A°-59A +10 , 
(a-2Ka 2) 7 (A-2}(A-4) (A -5) 
2(A—1-7680+12736A —7776A? + 20484 -94A*—63A°+10A°) , 
(A-2F(A-47 (A -5A -6) : 
2A(A-1)  _ 3(2016-2064A +3804? +1794? -67At +6A°) , 
(A =4Xa 6)” (A -3A -47 (A -6f (A -7) Í 
GA(A—1)-129024+136416A — 53408A? +9152A° —433At -59A° +6A°) ; 
(A -3A -4F(A-6F(A -7A -8) is 
2A(A-1)  51840-39456A+6312A°+1275A°-397A*+26À , 
(A -6}A -8) ? (A—5XA-—6)(A-8) (4-9) : 
2A(A -112856320 -14298624A + 62263684A? —1329216A° +135954A -361 1A% -4174$ +264) ; 
(A—4{A-5X{A-6) (A -87 (A-9YXA-—10) a 


n=0>b"(g)=(1-A} 


n=1b;(4)=(G-A) 


n=2b;(g)=(5-AŸ 


n=3— b"(q)=(7-AŸ 


En fait il s'agit des mêmes valeurs propres que pour les développements pairs mais en réalisant la 
substitution A->1- A. 


Solutions algébriques des équations différentielles de Mathieu associées 


A l'occasion du point précédent, nous avons vu que les développement des valeurs caractéristiques 
lorsque q proche de 0 se simplifiaient grandement pour la valeur A=-1/2 : 


A = N _ j 

asla) TE or 
2A-1) 1 

A(a)=(1+A) +2 

añ(g)=(1+ A) 4 +29 


Il s'agit des valeurs exactes des caractéristiques de l'équation de Mathieu associées liée à des 
solutions algébriques (voir le document sur l'équation des ondes sphéroïdales avec u de valeur 
=] 2 2 
Fu kan u= k q=% PEAN RYA E 
entière) : Bo 2 2 a A: . Les solutions de l'équation 
intermédiaire de Mathieu associée sont “®)=4 Cos2/aCos(s)}+ B Sinl2 a Cos(s) Je rappelle le lien avec 


entre les fonctions de Mathieu associées intermédiaires et les fonctions sphéroïdales : 


(il -)22@ 2x 40) Ç H 4 qû x?) A p6)=0 


P O EE E EEN A EE E E E O 


Insérons dans l'équation intermédiaire de Mathieu associée un « ansatz » de la forme: 
à > l=u-l1 

u()=e "5% p (cos) p=l-A P()= dax à,,=1 

2 1=0 , où Pu est un polynôme de degré u-1 en 

Cos(t). La construction des solutions est identique à celle des solutions algébriques l'équation des 


ondes sphéroïdales, à savoir, pour ,;,-2 A= Jen - +24, la solution est de la forme : 


u(t)= g?fa coslo) (cn | avec @ solution de P(@)=ø@(@-2)+16q=0 


4 Ja 


Pour y=3 A=- a=0++924; la Solution est de la forme : 
2 4 


4 Ja 32 q 


u(t) = et 0 (cwt) H- Cos(t) H dv) avec & solution de P(@)= © lo 2o 6) +64q lo 4)=0 


Pour 4-4 A=-7 a=œ+Żl+2q: la solution est de la forme : 
2 4 


(= 6 0 | os ia 2 argy O Ra oa ENEN oa 2e) 
4q 32 q 384 q/q 


avec © solution de P,(w)=«(o-2(o-6(o-12)+32q (5 œ -66 © +180)=0 


La récurrence générale que l'on peut établir sur ces solutions algébriques est identique à celle 
obtenue sur les fonctions sphéroïdales. Les valeurs propres w racines du polynôme s'obtiennent en 
remplaçant dans l'équation différentielle sphéroïdale un polynôme de la forme : a, +a, x+ a,x? 


pour ne retenir que le terme de degré 0 qui doit s'annuler également, soit la condition sur les trois 
coefficients : a(o+u(i-u))+2(qV2 a +a,)=o- Cette expression constitue un polynôme de 
degré m en À. En prenant le motif suivant intégré dans l'équation équivalente différentielle des 
ondes sphéroïdales : ax! +a, x"! +a, x"? +a, x +a, axt +a xs. Puis en recherchant le 


terme de puissance de x comportant 4 termes successifs, il y en a trois mais qui se déduisent l'un de 
l'autre par translation d'indice puis en annulant l'un des trois termes, on obtient la relation de 
récurrence suivante : 


=1 a = 


i œo 
u—2 4/q 
-4i Jq(l+1-u)a, -(2+P -0+1 (3-2u)-3u + u? )a,ı +(1+2Y4i Ja a: +(+3)a.)=0 
do = aÿ(@) polynôme de degré u-— 1l en © 


a =a (œ) polynôme de degré 1 —-2 en œ 
> aÿ(@ + ul — u))+ 2 (2i Va a + a,)= O polynôme de degré u en © 


Pour l'exemple les premiers termes du développement s'écrivent : 


u(t)= 62 fa 0 Cos"(t}+i Cos“?(r) H 2 G=o)+16(1-u)g Cos“*(r) +i 16qø 6 3u) t324 W 1) olo 2\o 6) Cost (t)+.... 
4Jq 32q 384 qq 
u>0 


Correspondance entre les valeurs caractéristiques des équations différentielles de Mathieu 
associées et des équations différentielles de Mathieu associées intermédiaires 


On a donc essentiellement 6 types d'équations différentielles reliées entre-elles, et pour chacun des 
deux groupes 2x4 solutions, notons les comme ceci : 


w''(t) +2A Cotan(t) w' (+ lag” (A,q)- 4° +2q Cos(2t)) w(t )= 0> ccei/\ ( t ,q)= Fa q) Cos( (21 t) 


w(t) +2A Cotan(t) w' (t)+ (as, (A, q)- A +2q Cos(21)) w(t) =0 > ccei,,., (t, q)= AA (q) Cos s((21 + 1) t) 


w"'(#)+2(1— A) Cotan(t) w(t)+ fac” (1 A, q)-(1- A} +24 Cos(21)) we) = 0 > ceci (t a)= D 45 (a) Cos(21 1) 
w (Q+ 201 À) Corali) w (+ asali Aa) 0 -AF +29 Cosa) le) =0 => ceii lesa)= F Ala) Coser +)1) 
Bo e(O) as,(a.g)-24 Coster) MAT (9 = 0 -> ceeste) (sn O) D Aila) Costa) 
sf as. (a9)-24 Coste) FET] eO =0 > cee alea) (S O S 43 (0) Coser) 
e(a (2s, a0)-24 Can) NA (= 0 > seh alia) (su (0) à È Bil) Siler +2)) 
e(a bsalasa)-24 Coster) MA Jel) = 0> exe (a) (su) E 83, O ser) 
0 
w(D=2A Tan) de 2A(a) Sin(2r+ 09 
w'(t)-20-A) Tanl) w (t)+ (a3;"(1 A,4)-(1- A -2qCos(21))n() = 0 seek" (e, SA a) Cos(21 1) 
(201 A) Tanl) w (o) hzl- A a) AT -24Cos21)) wli) =0 sera) D BA (a)Sin((21+ 09 
D T A a gl) = 0 => sees, (q) = (cos 9} X 4l) Cosl2r 9) 
g"(e)+| aza lA q)+ 24 Cos(2r) ael- 0 sce, (t,q)= (Cos? (e)) 2 FAO) )Cos((21 +1)1) 
eO (r ag) 24 Can) (= 0> sf, 2629) (Cos (0)? 83.4) Sr +2)) 
eO (bilas) 29 Coster) MA 9-0 se (a) (Co È BA la) Sir) 


Par ailleurs notons bien que pour les équations différentielles de Mathieu associées intermédiaires, 
je n'ai réussi à faire converger la récurrence à 4 termes que pour les valeurs de u demi-entières, 
donc pour des valeurs de A entières. 


Les formules de liaison de ces diverses valeurs caractéristiques sont les suivantes : 


dy (A, q) = a (A.-q) 
azn (A, q) = Aypa (A.-q) 


x= Cos(t) j 
Bu (A, q) = a (1 = A.,-q) 
bis (A, q) = CPE (1 z A,-q) 
a3, (A, q)= ax" (A-4) 

x= Sin(t) Arn (A, q)= dr (I-A.-q) 


bi,s2 (A, q) = by, (1 z A,-q) 
bou (A, q) = bya (A,-q) 


Construction de solution de l'équation de Mathieu associée modifiée 


Cette construction est très liée à celle des fonctions d'onde sphéroïdales radiales par un 
développement sous forme de fonctions de Bessel sphériques : 


H 
xell,+0]  S#(x,72)= (i = =) Srt Vaal x) i=12,3,4 
X k=—0 


Fonctions de Bessel Sphériques SÈ (x) fË = j, y,h®, hO 


Cette construction conduit à une récurrence similaire à celles des fonctions angulaires sous cette 
forme: 


Arfa (7° )+ (z, art, (7° )+ Ces (7°)= 0 ke {— 00,...,—1,...,0,1...,+00} 
Le 2(v+u+2k—1)v+u+2k) 

i (2v +4k -32v +4k —1) 
o (V +2k\v+2k+1)+ u° -1 

(2v +4k —1)(2v + 4k +3) 
x 2(V-u+2k+1{v-u+2k+2) 
(2v +4k +32v +4k +5) 

1 1 1 1 

2v +4k -3° 2v +4k—1° 2v +4k+3°2v+4k+5 


B, =(v+2k\v +2k+1)-2y 


Pôles possibles sur 


Cette récurrence présente également des pôles pour les valeurs de v demi-entières. Aussi il est 
également nécessaire d'envisager une autre forme de construction des fonctions sphéroïdales. 


Exprimons par exemple la fonction sphéroïdale radiale de première espèce : 


D'une manière tout aussi similaire à la construction des solutions radiales de l'équation de Mathieu 
modifiées, il est envisageable d'introduire l'équation différentielle des fonctions associées de 
Mathieu modifiées. Pour cela il suffit d'introduire ce que j'appellerais la forme de Liouville 
« hyperbolique » de l'équation différentielle sphéroïdale. Elle est obtenue avec le changement de 
variable x=Cosh(t) : 


2 
2 
1—x 


q-e) Ee aS orh x?) 


po =0 x= Cosh(t) 
y(x) = (Sink (E)E w(0 = (Cosh? (t)-1) = we) 


> 4 A = 5+ u q = — > w''(t)+ 2A Cotanh(t) w' (t)- (a — N —2q Cosh(2t)) w(t) =0 


2 2 
Y 2 y 1 
——U— =D+—+— 
SAT 2 ‘4 


a=@+ A + 
Une forme équivalente est obtenue avec y=i Sinh(t) : 


(i a) E ay Bia | Ç wis L }«-0 x = i Sinh(t) 


(x) = (Cosh? (£) wle) = (Sinh? (e)+ 1 = we) 
= ASE + q=—-— — w" (t)+2A Tanh(t)w'(t)- la- À +2q Cosh(2t))w(t)=0 
Z d 


a=0+% +- 
2 4 


En transcrivant le procédé de l'article de 1923 de E.L.Ince, « Associated Mathieu Functions », le 
point de départ est la forme de Liouville « hyperbolique » suivante : 


= Cosh(t) 
ya)= (Sin (9) e) => gl) = (Cosh? l) -1 wlr) E 
ie e + g QE -24 Cosnlzi)+ MA: Peo- o 
a=@+—+ 
4 2 
Une forme équivalente est obtenue avec y=i Sinh(t) : 
X=i Sinh(t) 
» = (Cosh? (0) g)> g(t) = (Sinh?) +1) wr) i 
ee > g" ()-(a-24 Cosni)- AE 20- o 
SERKA 
a ES 


Si l'on doit se référer aux équations différentielles de Mathieu associées, avec leur forme 
intermédiaires : 


xr 


w"'(t +2A Cotan(t) w' (t)+ (a — N —-2q Cos(2r)) w(t) =0 
w''(t)— 2A Tan(t) w'(t)+ (a — N +2q Cos(2r)) w(t) =0 


g"'(e)+ (a D AG W) Sme 


Sin? (t) 


g''(e)+ (a —2q Cos(2s)- AA ©) un 


Cos’ (t) 


On voit qu'il suffit d'effectuer la transformation t->i t, pour obtenir les équation différentielles de 
Mathieu associées modifiées. 


Aussi une telle construction est donc en théorie possible avec toutes les formes de développement 
en sinus hyperbolique, par exemple : 


x= Cos(t) x= Cosh(t) 
g''()+ Ç — 2q Cos(2t)— ay =0 g' '(z)— Ç — 2q Cosh(2t)+ MAR) =0 
ge) = (Sin (OF = A,,Cos(21 t) ge) = (Sin (OF = A,,Cosh(21 t) 


gle) = (Sin (OF > AnCos((2i+1)r) 7 |e()= (Sinn (0) > As Cosh((21+1)£) 
ge) = (Sin (O? > BnSin((2i +1)r) g(t) = (Sink? (OJ? 2 B, Sinh((21+1)t) 


gle) = (Sin? (O? pa Ba Sin((21+2)t)  |g()=(Sinn?(r))* > B.,.Sinh((21 +2)r) 


x= Sin(t) x=i Sinh(t) 
g''(t)+ (a — 2q Cos(2t) aay] g(t) =0 g'(t)— (a — 2q Cosh(2t) 2 g(t) =0 
ge) = (CoP (OF > A,,Cos(21 t) ge) = (Cos (OF = A,,Cosh(21 t) 


— 


g) = (Cos A? D Aru Cos (2 +1)z) 


get) = (Cos O° D Ar Cosh(@1 +1)z) 


gle) = (Cos (F _ Ba Sin((21+1)£) gle) = (Cos (OF 2 Ba Sinh((21+1)£) 
gle) = (Cos (O)? 2 Bao Sin((21+2)t) gle) = (Cosh (N? > B, Sinh((21+2)t) 


Le gros inconvénient de ces constructions c'est qu'elles ne convergent pas rapidement, voir pas du 
tout au delà d'un certain rayon de convergence. 


La question que l'on se pose immédiatement : existe-t-il une autre forme de construction pour les 
fonctions de Mathieu associées modifiées dont la convergence serait plus rapide. Je pense 


notamment à des constructions comme celles des fonctions de Mathieu modifiées à l'aide de 
fonctions de Bessel. 


Commençons par l'équation "'"(t)+2A Cotanh(t) w'(t)— (a— À -2q Cosh(2t))w{t)=0 et écrivons sa 


version algébrique à l'aide du changement de variable £=24q Cosh(t ) il vient : 
(£2-4q)w"(£)+(1+2A)E w'(8)+ (87 + À -a-2q)w(£)=0 
La forme « asymptotique de cette «équation différentielle est la suivante : 
En" (E)+(1+2A)E w(E)+(E7 + A —a-29)n(E)= 0 


Avec le changement de fonction : m(E)= EVE), il vient : 


EO OG a—2q aq NN} ue) 0 


Cette équation a pour solution asymptotique en € grand : 


E? v'(e)+8 v(E)+ (65 -a-2q)E)=0 > E)= J (€) 
On peut donc proposer de construire une solution de l'équation : 
w''(1)+2A Cotanh(r) w'(t)— (a — N -2q Cosh(2r)) w()= 0 


Sous la forme d'une série : (= S 7Y)  %()=(Cosa())"J (2/4 Cosh) 


Comme les fonctions Y(£) sont solutions de l'équation différentielle : 
E y (E)+(1+2)8 y (E)E + A7 -7)r,(6)=0 
Et que la forme algébrique de l'équation différentielle des fonctions w est la suivante : 
Ew,ë}= (6? —-4g)w"(e)+(1+2A)€ w(E)+ (E? +A -a-29)w(£)= 0 
Il vient : 


E{w,£}=0= Dre —4g)r"(E)+(1+2A)€6 7'(E)+(E2+ A7 -a -2q)Y,(£))=0 


Or E? Yy (E)= -0+ 24) PE- E+N -PES Dr -a-2a)rlE)= 45r" E) 


Or l'expression * (E) peut être « linéarisée» à l'aide de ses proches voisins : Y- (ELY (EY. (E ) 
dont voici l'expression : 


A= aU FX a 


B, (CS AJl SCE) > Y (E)=C, Yna(E)+ BY (E)+ 4 YCE) 


C, taa 14 2 
4 I +1 


Il vient donc finalement la loi de récurrence à trois termes lorsque l'on annule chaque terme 16) 
dans la sommation : 


((—-2+A)7-1+A) ~en) À U+2-AX+1-A), zô 


H P2 
Gene Ce (a Pee G+2X1+1) 


Cette récurrence redonne bien la récurrence du développement des fonctions de Mathieu modifiées 


TE > yJ, (2/4 Cosh) 


lorsque A=0, de la forme : , qui est celle des fonctions de Mathieu (à un 


.2 AN 
changement de signe près) : 1712* (a-p )ri +a Yn = o: 


Il me semble dans ces conditions que nous nageons dans le bonheur, mais tout s'effondre en réalité 
car le premier terme du développement hélas ne respecte pas la règle de « linéarisation » : en effet 
cette expressions n'est pas définie pour 1=0 et 1-1. Cela se concrétise également par l'apparition 
dans la récurrence à trois termes de pôles justement pour l=-1, 1=1,1=-2 et 1=2. (Une récurrence de 
forme similaire est étudiée dans l'article de 1950 de R. Campbell « CONTRIBUTION A L'ÉTUDE DES 
SOLUTIONS DE L'ÉQUATION DE MATHIEU ASSOCIÉE », pages 214 et suivantes, mais elle s'avère 
justement fausse de part une petite erreur de calcul, mais aussi pour avoir négligé de tester 
l'application de cette récurrence au premiers termes du développement. 


Reprenons ainsi le premier terme du développement 1-0 par exemple : 


WEOE EAE) =ne E or 10 EE) 


ua, A(A +1) 


P(E)= EE) 17," E) = -Enn (E)+2AE die J (E)= 7, A 2A) z 


1— 


2 2A y (E) 1+2A 


A(A +1) 
= 


»,(£)+ nhé) 


En injectant ce premier terme dans l'équation différentielle, il vient : 


(-a-29)Y,(E)= 447," (E)= aar -= Arr AU br 1()) 


5? 


A(A + ns 
a 


S -yola +44) Y (E)+ 240 +24) (6) = 7,49 — (E) 


Si maintenant on introduit l'équation inhomogène suivante : 


ASE) = 74005 (0) 
ZOR Erne) WE) 


Elw, E} = (E? -4g)#"'(E)+(1+2A)€ F'(E)+(E2 + À a -2q)HE)=-7,4q ET A(A +1) 


mr nE) 


> Efi e]=0 (6 44) P(E) +0 2A)8 PEHE + A a 29) (€) 0 


Est-il possible de construire un développement effectif de la solution de l'équation différentielle de 
Mathieu associée modifiée intermédiaire ? Pour cela pouvons-nous construire simplement une 
solution particulière de l'équation inhomogène : 


ag AA +1) +1) 


a M) 


(E? 4q) W (E)+(1+2A) E WE) E + A -a-2q) WE) = -y 


A ce stade je n'ai pas trouvé de solution particulière à cette équation inhomogène. 


Y(E) 


Le fait d'aboutir à une impasse pour la construction des solutions de la forme 
w(t)= aT D È, J2 Ja Cosh(t)) pour les solutions périodiques est somme toute normale, car 


ce sont en réalité les solutions radiales de l'équation des ondes sphéroïdales. Pour cela reprenons 
une ra non nécessairement périodique pour l'équation de Mathieu associée modifiée : 


WP, 1)= ST Va Coho) 
Cos TT) — , elle suit la récurrence suivante : 
(5 +21-2+AX5+2/-1+A) z 2 A(A—1) (5 +2/+2-AX5+27+1—A) 
= = m+|a @ H 21) 2q -= z 1 1+1 
© +2127 +21—1) (5 +27) —1 (5° +21+2Y7 +21+1) 
2 72 
: Ke ty q= aps gamer pe , 

Selon la valeur des paramètres : 2 4 2 2 4, cette récurrence 


s'écrit : 4Y + (B, -@)y, + Cri = 0: avec les valeurs des coefficients : 


2 v+ u+21—1)\(v + u+21) 


À, = 
: (2v +41—3)\(2v +41-1) 
j i 2: 
=(v+2/)v+1+2/)-2y AU) CE ET ae 
(2v +41-1Y2v +41 +3) 
a > Vout 1Xv-—u+21+2) 
d (2v +41 +32v +41 +5) 


C'est bien la même récurrence obtenue pour la construction des fonctions d'ondes sphéroïdales 
radiales, et de plus le développement s'écrit effectivement de manière similaire. 


x = Cosh(t) 7 =2V/q =v 


1 1=+00 


(Cosh? (t)- 1)? sO” (x,y?) > ri Jo Jq Cosh(t)}< w(t) 


Cosh^ (t) A v,k® V+21 


HO erere Ho" Ja Cosh(r)) 


Il n'y a donc rien d'étonnant à ce que la récurrence ne marche pas pour les cas Y <tier puisque 
pour les fonctions sphéroïdales radiales la récurrence n'est pas définie pour v demi-entier. 


Conclusion : on peut effectivement construire des solutions non periodiques de l'équation de 
. s rs PIENE wẸ.t)= Co: TÈS Iosia Cosh) | 
Mathieu associée modifiée intermédiaire de la forme 1) io 
la seule condition que Ÿ ne soit pas un entier. Mais cela n'est pas vraiment un duia 
supplémentaire puisque il s'agit de fonctions apparentées aux fonctions sphéroïdales radiales. 


Par contre lorsque Ÿ est entier, ce qui correspond à des solutions périodiques pour v demi- 
1 l=+0 
we) = 2 Sy, 2a Coshle)) 
entier, la construction de sous forme de ce développement Cosh^ (t) 
n'est pas possible. 


Tableau synoptique de la construction de solutions de première espèce de l'équation de Mathieu 
associée modifiée et de l'équation de Mathieu associée modifiée intermédiaire à l'aide des 


fonctions de Bessel sphériques 


La récurrence générale des coefficients 


des développements des fonctions d'ondes sphéroïdales 


Abi r2 + Bbi ra T Cbi k + D, be + Ebi k2 = 0 
A = y (0° 1) (~v-u+k-I\v-u+k) 
(2v +2k -32v +24 —1) 
EEP „etiv = u+k) 


C, =(v + kv +k +1)+27° la? 


2v+2k-1 


(v+kv+k+1)+ u’ -1 
(2v +24 —1)(2v +2k +3) 


1) 


D -2a (v+k+1\v+u+k+1) 
$ 4 2v +2k+3 

iare a 
Š (2v +2k +32v +2k +5) 


La récurrence générale des coefficients des développements des fonctions de Mathieu associées et 
des fonctions de Mathieu associées intermédiaires : 


A, bôi +B, bår +C, bô, +D, bôi +E, bô 50  y=24q a=4+2q+2 u=}-A 
E CR ee do me 
a Cv +2k-3X2v +2k—1) 
(2v +2k)\(2v +2k—-1+2A) 
B == 
emaa 2v+2k—1 
C (vers) +2 a? +8 qla? —1) A(A —1) 3 
LE 2 z 1 (2v +2k -12v +2k +3) 
(2v +2k +2 2v +2k+3-24A) 
D = 
waana 2v +2k+3 
E (a Tea aE E 
A (v +2k +3 2v +2k+5) 


Domaine 


AS Cosh(t) 


Développement des solutions radiales de 


l'équation des ondes sphéroïdales 


herezio art-e)- a 


Développement des solutions de l'équation de Mathieu associée modifiée 


e"O-[a —2 q Cosh(2t)+ D Le) Oo et lien des 


solutions et des paramètres : 


y=2/a A : ua 2+7+24 g) = (Sin? (NF (€) 


Développement des solutions de l'équation de Mathieu associée 
modifiée intermédiaire : 

w"(t)+2A Cotanh(t)w'(t)- (a -A -2q Cosh(2t))w(t)=0 
et lien des solutions et des paramètres : 


y 2/4 A u a A+2+29 we) = (Sin (9)? 80) 


2 


Récurrence et 
paramètre q 


E= Cosh(t)>1 


yl) = A Pau (0) sal NE) 


g(t) = (Sint (0) De (0) Jak Va Sinh(t)) 


1 A k=to 
w(t) — (Sink? ()} 2 + DÉaPhiox (0) Jv+2k Va Sinh(t)) 
k= 


Récurrence à trois 
termes et a=0 


é= Cosh(t)>1 >l 


k=-% 
D Sp, P 2k PA Ohaba] 
k=-0 


k=+00 


(= (Sinh (9) +? Coshle) SH, Paa (O) jaa Sinne) 


k=- 


Récurrence à trois 
termes et a=0 


ES Cosh(t)>1 


OERS | ritt- a) 


gli) = COH (SOE Terra (O) Seala Sne) 
[EE siet OF Teri Va Cosh(t) _ a) 


v(t )= (Cosh(t)- o Jaa ( Cosh(t) -a)) 


Récurrence à cinq 
termes et a 
quelconque 


E = Cosh(t)>1 


eA 
olti 7 


lesa} J £ I-v-k-afll+v+k-y 


jil 11) 


gt ÿ | "a jaloa (Comi) 


Sink’() i T(-v-k-uE(l+v+k- u) 


1, k= b! 


w{t)= (Cosh) +a) lsi (0) L l n Ja ba (Cosh(}+a) 


A T(-v-k-uEll+v+k- y) 


Récurrence à cinq 
termes et a 
quelconque 


E = Cosh(t)>1 


SR 

1 2 +0 | 

l1- E: ) Do ne (y é) 
kao 


g(t)= (Tanh? OE (Sinh? (} De in Va Coshlt 


w(t) = (Cosh) ^ Se Jra (Va Cosh(r) 


Récurrence à trois 
termes et a=0 


E = Cosh(t)>1 


1\2+ br, 
-b-k A T(-v-2k-u(1+v+2k- u) 


aly d 800= a(o} 


1A 1 k=+0 H 
go z b, 


2 (Sin (o)} > 


I(-v-2k-u(1+v+2k- u 


ji (Ja cosita) 


Mo) =(Cos) - ul 5 ba Pa (Costo) 


Récurrence à trois 
termes et a=0 


Cosk(t)) Æ T(-v-2k-uE(l+v+2k- u) 
éz Cosh(t ) >l b a Récurrence à trois 
o To Daa Jal y(é=+1)) gl [eta] ` (Sink?( 0): > be ja Va ( Cosh(t )+1) w(t) = (Cosh(t)+ DT CT Cosh(t)+1 ) termes et a=-1 
é= Cosh(t) >l Ee Hi 75 = (Cosh 5 2 Ja (Coshlt Récurrence à trois 
Efsun EE wo Bebia OA Benbge EE 


Domaine Développement des solutions radiales de | Développement des solutions de l'équation de Mathieu associée modifiée | Développement des solutions de l'équation de Mathieu associée | Récurrence et 
l'équation des ondes sphéroïdales ( ) es e ) A( A— 1) ( ) Oi ret lien: ides modifiée intermédiaire : paramètre q 
2 "(z)— a A t t)= 
E = Cosh(t) 20-260 arhei- g 052a COSA Sinn e) E w" (t)+2A Cotanh(t) w'(t)- (a-~ -2q Cosh(2t))w(t)=0 
á i í 1-€ POS solutions et des paramètres : et lien des solutions et des paramètres : 
7=2/q A : u Oa=à+ " +29 gle)= (Sin (OF y) y=2ļJq A=- a=1+ f +2q w(e)=(Sink (e)? g4) 
7. gl e b“ 75 1 Ex 4 Lee br, Récurrence à trois 
@{ ) E 6) | a (GLS ok bi u(t)=(Coshlt)-1} r jala (Coshlt)+1)) | # tà=-1 
d ë+1 È TCV kalv rku) gli)= Cos =1 (Sinh o} 2 p jt (Ca) Æ TCv-k-u(l+v+k-4) 1 | ermes et q 
€ >1 u ae i 1, k b! Récurrence à trois 
+l bsx 7 z Cosh(t)+1 £2 LE bii ; t}=(Coshlr)+1}7* vá java (Cosh(t)-1 t a=1 
"E- . 3 > Rp CENTRAL 1) eo- (Sime (} D a a (fr (Cast) 0 ©) ) À T(-v-k-uT(1+v+4-u) | va o) ) CE 
Ez Cosh(t)>1 ) >l i Récurrence à trois 


=> bii > (+ Epal E) 


u DA 1 D 
CURE =) Jv+k (Jq Cosh(t) 


w(r)= (Sinh? Oa à iS mr (+ 


k=—o 


p) ja Va Cosh(s) 


Cosh\t 


termes et a=1 ou 
a=-1 


Dans le cas où A=0 (fonction de Mathieu) alors la récurrence générale s'écrit : 


Ar bise 2 +B, be x 1+C% br: +D, b, p1 +E; bjr: =0 


A, =q læ? -1) 
B, = à q (2v +2k) 


2 
C, (ver) +2qa° — 


, = a q (2v+2k+2) 
E, =q (a° -1) 


Lorsque a=0 alors la récurrence saute de deux en deux sur les indices k : 


i 1Ÿ 
Gus =(- 1) bo, 2 > 4 a, a2 H (e 2k + >) a] a or +a D? araa =0 


1y 
RE =(- 1) bo, au 7 4 a, a H (e t2k+1+ >) a] CTA +q a =0 


Cela correspond donc bien à la récurrence des coefficients du développement de la fonction de 
Mathieu angulaire à un changement de paramètre près. 


Le tableau synoptique suivant reprend les développements possibles de la fonction de Mathieu 
radiale (ou modifiée) : 


Développement des solutions de l'équation de Mathieu modifiée | Récurrence et paramètre a 
g"'(t)- (a —2q Cosh(2r))g(r) = 0 


k=+00 1 
er) È BaPZal0)I Va Sinne) 
k=—%0 2 


1 


g(t) = Cotanh(t OPA by, Par (0) Jys oa Sinh(t) 


k=- 


Récurrence à trois termes et a=0 


Récurrence à trois termes et a=0 


E) Récurrence à cinq termes et à quelconque 


gl) = or (Va (Cosn(s 
t3 


Récurrence à cinq termes et à quelconque 


el)= Gt S: bs 3,04 (Coshle)+a) 


Coshlt)+a Æ T(-v-k-uFll+v+k-u) 


k=+00 
g()= X Bag (Va CosA(r) 
k=—0 2 


Récurrence à trois termes et a=0 


Récurrence à trois termes et a=0 


g(t)= Tanh(r) > b? a(3 sv Jr (Va Cosh(t)) 


Récurrence à trois termes et a=-1 


el) = È B, (Ja (Cosh(s)+1) 


) Récurrence à trois termes et a=1 


gle) = È be (Va (Cosh(s)- 


er Récurrence à trois termes et a=-1 


$ to 


Ye E svk} J, CVa (Coshle)+ 1) 


Récurrence à trois termes et a=1 


g0)= e T Èe [jevreja (Cosh(t)-1)) 


Cosh ()- 


1 Récurrence à trois termes et a=1 ou a=-1 
a0) = (Tan OS mr +]... Va Conte) 


k=—0 


Pour lesquels on reconnaît déjà deux développements connus des fonctions radiales : 


k=+00 
g()= X beag (Va CosA(r) 
k=—co 


+—+24k 
2 


g(t)= Tanh(t) Seal svp, a Ra Cos(r) 


Or les expressions des fonctions de Ferrer sont : 
i Cos (v + >) AreCos(a)) 
1 2 2 
(= 2 | 
4 (x) 


so eme) (fret) eme) 


2(1-x2} (i-x2}s 


Cela conduit à évaluer les deux expressions comme suit : 


1 
P2(0)= [2 Cos[ É 
TT TT 2 
ArcCos(0) = — 


T ho Ee Ee) 


Pa (0)= 2 Conf ? É | à) | tx }-( 1) P2(0) 
VOREG | aket) sedye) | een o 


NI =e 
NL 7 
KL 


Soit : 


z k=+0 12 . : 
Le développement (= 5° 5°. P2 (0) A e Ja Sinh(t)) s'écrit donc : 
k=-00 Vegte 


(2 Ja Sinhle))= a la CII, ra eVa Sinh(£)) 


v+— 724 


s()=P2(0) s Baly 


et le développement o(t)= Cotanh(t jE b? LP: (0 Fe k Va Sinh(e ) s'écrit sous la forme : 


1 


0-20 ) Cotanh(t D po, y +542) J, (Va Sinh(o))o Cotanh(t p3 b (- yf +342) J, Na Sinh(o) 
pel 
2 


On reconnaît dans ces deux développements clairement des développements connus des fonctions 
de Mathieu modifiées. 


Développement de la fonction d'onde sphéroïdale radiale valable pour les valeurs u demi-entières 
et des fonctions de Mathieu associées modifiées et intermédiaire valables pour A entier 


Au paragraphe 3.65 de l'ouvrage «Sphaeroid funktionen » de J.Meixner et FW.Schafke en formule 
(49) page 293, les auteurs proposent un développement des deux solutions radiales de troisième et 
quatrième espèce qui est valable quelque soit les valeurs de u, pour z dans tous le plan complexe 
hormis le demi-axe des réels de -c à 1. Il est donc valable pour les valeurs réelles de z de 1 à +% : 


En injectant cette forme dans l'équation des ondes sphéroïdales et si la forme proposée doit 
respecter la condition de normalisation des fonctions radiales à l'infini, alors les coefficients du 
développement suivent une relation de récurrence à 4 termes : 


f y’ (l-u\l-u-1)4, +47- u)A" a +i 0+- A + (+04, 0 
Æ =4,=0 Æo =l 


Dans le cas où u est demi-entier alors cette récurrence est parfaitement définie pour toute valeurs 
de l'indice |, et la forme des fonctions radiales de première et deuxième espèce définie comme suit 
est donc parfaitement définie pour pu demi-entier et v demi-entier, cas qui nous occupe 
particulièrement car jusqu'à présent aucun développement n'était défini pour ce cas. 


2 e : Z (2) 2 5 5) 2 e i Z 
LS Ea| P È z 2% Car 


Pour les fonctions de Mathieu associées modifiées respectant l'équation différentielle : 


e"0-a-2 q Cosh(2t)+ saho- 0 


Sinh’ (t) 
Avec comme lien des solutions et des paramètres avec l'équation différentielle des ondes 
sphéroïdales : 
1 
png Agen a=à+i+2g glo=(Sm (0) le) E= Cost) 


La récurrence s'écrira comme suit : 


4q{(21+2A-1f21+2A-3)4", , +8 q (21+2A-1)A^ + L (2+1)+ a+ 2q- a +(1+1)4},, =0 


À, a AS =0 A =1 
ifada Cosi) z) 


A l l=+0 zi 
gilt) = Re (Sink? (£) Cosh? (É Cosh(t) 


) 2./q Cosh(t) 24 > 4i Ja) 


` 1. ifaya Cols) a F 
8.021 (smel iO | Co) 


2,/q Cosh(t) 24 Caia] 


Pour les fonctions de Mathieu associées modifiées intermédiaires respectant l'équation 
différentielle : 


w" (t)+2A Cotanh(t) w'(t)- (a -X -2q Cosh(2t)) w(t)= 0 


Avec comme lien des solutions et des paramètres avec l'équation différentielle de Mathieu 
associée : 


2 
2 


u(t)= (sinto) gli) 


(l-v -1)r 
2 


s E a OED Cos (zya Cosh(t)+ 
y- g a” (aJa) Cosh'(i) 


ig (1 =V ir 
FE (9) pe Sin Er Cosh(}+ == 
W; (e) z 


Tu 4 (Ja) Cosh'(s) 


Le cas des fonctions de Mathieu modifiée est obtenu pour A=0, soit : 


w'(t)-(a-2q Cosh(2t)) w(r)= 0 


Avec comme lien des solutions et des paramètres avec l'équation différentielle de Mathieu 
associée : 


4q (2/ -1)21- 352 +8q (21 = 14: «| (+ 1)+ ; +2q- a+ + (1 +1)4? = 
4, 4 = 0 4 =] 
Are À > fs (4 Ja) Cosh((t) 
z Sin U Cosh(t) fr) 
Er nn À 2° (4 q) Cosh'( (1) 


Cos (zya Cosh(t)+ eok) 


w(r)= 


